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information  storage  time  of  more  than  two  months,  and  a  framing  rate  of  40  Hz 
with  full  modulation  depth  (200  Hz  with  20%  modulation  depth).  The  theoretical 
details  of  operation  of  the  device  are  described  in  the  Applied  Optics  paper 
which  constitutes  Appendix  I  of  this  report.  Recent  advances  in  crystal  polish¬ 
ing  techniques  have  resulted  in  the  availability  of  thin,  high  quality  LiNb03 
crystals,  and  in  a  concommitant  increase  in  the  spatial  resolution  and  the 
optical  quality  of  the  device.  Spatial  resolutions  of  better  than  2  cycles/mm 
at  50%  contrast  (^10  cycles/mm  at  10%  contrast)  have  been  achieved.  Additionally, 
we  have  shown  theoretically  and  experimentally  how  the  MSLM,  functioning  as  a 
phase-only  modulator,  can  be  successfully  integrated  into  high  resolution  all 
optical  adaptive  systems  to  perform  wavefront  phase  compensation  for  applications 
such  as  atmospheric  optical  communication  and  imaging  through  randomly  fluctua¬ 
ting  media. 
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PREFACE 


Development  and  Applications  of  the  MicroChannel  Spatial  Light  Modulator 


This  final  report  summarizes  the  MIT  program  on  the  development  of 
the  optically-addressed  MicroChannel  Spatial  Light  Modulator  (MSLM).  The 
work  was  sponsored  by  the  Air  Force  Office  of  Scientific  Research  under 
Grant  No.  AFOSR-77-3328. 

Since  the  basic  theory  of  the  operation  of  the  MSLM  has  already  been 
published  in  open  literature,  only  a  brief  summary  is  given  in  the  body 
of  the  report.  The  reader  is  referred  to  Appendix  I  which  is  a  reproduc¬ 
tion  of  the  literature  on  the  theory  of  operation. 

'^>Vhe  fundamental  operating  characteristics  and  materials  limitations 
of  the  device  are  discussed.  The  role  of  secondary  electron  emission  in 
the  operation  of  the  device  is  stressed  and  some  of  the  write,  cycling, 
and  readout  modes  and  their  limitations  are  described.  In  addition,  the 
limitations  of  the  inherent  space-domain  image-processing  operations  of 
the  device  are  discussed.  In  the  final  section,  the  report  summarizes 
some  of  the  key  results  on  the  application  of  the  MSLM  to  the  problem  of 
high  resolution,  all-optical,  adaptive  wavefront  phase  compensation  for 
communication  and  imaging  through  the  earth's  atmosphere. 
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INTRODUCTION 

The  optically-addressed  microchannel  spatial  light  modulator 
(MSLM)is  a  versatile,  real-time  optical  signal-  and  image-processing 
device  that  exhibits  high  optical  sensitivity  and  high  framing  speed.  It 
consists  of  a  photocathode  and  a  microchannel  plate  (MCP)1  in  proximity  focus 
with  an  electro-optic  crystal  platef/^as^showrTTrr Fig.  1.  A  planar  electron 
acceleration  grid  is  often  placed  between  the  MCP  and  the  electro-optic  plate. 
The  electro-optic  plate  carries  a  high-resistivity  dielectric  mirror  on  one 
side  and  a  transparent  conducting  electrode  on  the  other.  . . 7  / 

During  operation  of  the  MSLM,  write  light  (coherent  or  incoherent) 
incident  on  the  photocathode  creates  an  electron  image  which  is  amplified 
by  the  MCP  and  proximity  focused  onto  the  surface  of  the  dielectric  mirror. 
The  resulting  surface  charge  distribution  creates  a  spatially  varying 
electric  field  that  modulates  the  refractive  index  of  the  electrooptic 


plate. 


Thus  the  readout  light,  which  makes  a  double  pass  through  the 


crystal,  is  spatially  phase  modulated.  To  erase  the  image, the  surface 
charge  distribution  must  be  erased.  A  variety  of  mechanisms  for  doing  so 
are  available,  depending  on  the  nature  of  the  charge  distribution  and  the 
optical,  transport  and  secondary  electron  emission  properties  of  the  surface 
in  which  the  charge  distribution  is  stored. 

The  versatility  of  the  MSLM  arises  from  the  architecture  of  the  device 
in viiich  the  functions  of  charge  generation,  multiplication,  transfer,  and 
storage,  as  well  as  that  of  light  modulation  are  performed  by  separate 
device  components.  This  design  permits  independent  control  of  the  above 
processes,  allows  the  device  to  be  tailored  for  specific  applications,  and 
leads  to  a  wide  variety  of  image  processing  operations  that  can  be  achieved 
by  operating  the  device  in  its  internal,  space-domain,  image  processing  mode. 


The  Optically  Addressed  MicroChannel  Spatial  Light  Modulator 


B 


» 
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In  general,  among  other  applications,  spatial  light  modulators  can 
be  used  as  phase-only  modulators,  amplitude-only  modulators,  intensity 
modulators,  Fourier-plane  spatial  filters,  input  plane  transducers,  and 
space-domain  processors.  Through  proper  choice  of  materials  and  confi¬ 
gurations,  the  MSLM  can  be  tailored  to  accomplish  any  of  these  functions. 
This  report  (1)  describes  the  basic  operating  modes  and  principles  of  the 
MSLM,  (2)  illustrates  how  these  relate  to  some  of  the  important  image 
processing  operations  that  can  be  achieved  when  the  device  is  operated  as 
a  space-domain  processor,  (3)  discusses  how  the  various  materials  properties 
impact  on  and  limit  device  performance,  and  (4)  explores  the  application  of 
the  MSLM  in  an  all-optical  closed-loop,  real-time  adaptive  spatial  phase 
compensation  system  that  could  prove  useful  for  communications  and  imaging 
through  scattering  media.  Earlier  stages  of  development  of  the  MSLM  have 
been  reported  in  references  2  and  3  (Appendices  I  and  II). 

OPERATIONAL  DETAILS 

Figure  2  is  a  sketch  showing  the  dimensions  and  voltages  of  the  device 
components.  Under  typical  biasing  conditions,  the  MCP  output  electrode  is 
held  at  ground  potential,  the  grid  at  Va  (~2kV)and  the  crystal  electrode  at 
a  potential  Vb  which  is  varied  during  operation.  The  voltage  across  the 
gap  between  the  grid  and  the  crystal  surface  is  Vg  and  the  voltage  across 
the  crystal  is  Vx-  Thus, 

\  ■  va  +  vg  +  vx  • 

and  because  the  capacitance  of  the  crystal  is  generally  much  larger  than 
that  of  the  gap,  the  changes  that  are  impressed  on  V5  during  operation  are 
coupled  primarily  across  the  gap. 
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Secondary  Electron  Emission 

In  some  MSLMs,  the  dynamics  of  the  manipulation  of  the  surface  charge 
distribution  are  determined  by  the  transport  and  secondary  emission  properties 
of  the  surface  in  which  the  charge  distribution  is  stored.  In  these  devices, 
it  is  important  to  note  that  microchannel  plates  are  not  monoenergetic  sources 
of  electrons.  The  output  electron  current  from  a  typical  MCP  has  a  continuous 


energy  distribution  n  (V.,)  that  can  range  between  zero  and  eV  ,  where  -V 

V  l\  III  Ml 

(Vm  >  0)  is  the  MCP  bias  voltage,  and  e  is  the  electronic  charge.  The 
distribution  is  weighted  heavily  in  favor  of  the  low  energy  (0  to  >50eV)  elec¬ 


trons.4’  5  Here,  nQ( V^)dVj<,  is  the  number  of  MCP  output  electrons  per  unit 
area  per  unit  time  with  energy  between  eV^  and  e(Vk  +  dVk).  The  MCP  out¬ 
put  current  density  JQ  can  therefore  be  written  as 


m 

Jo  -  (  e"o(IwVdVk  C2) 

0 

In  the  linear  operating  regions  of  the  MCP  and  the  photocathode, 
nn(I  ,V. )  scales  linearly  with  the  write  light  intensity  I  .  In  terms  of 

O  W  lx  Vi/ 

the  actual  device  parameters,  JQ  can  be  obtained  from  the  expression 


Jo  =  IwneG/hv’  (3) 

where  n  is  the  quantum  efficiency  of  the  photocathode,  G  is  the  gain 
of  the  MCP,  h  is  Planck's  constant,  and  v  is  the  frequency  of  the  write  light. 
The  primary  current  Jp  is  defined  at  the  surface  of  the  crystal ,  and  in  a 
gridless  device  Jp  =  when  none  of  the  MCP  electrons  is  locked  out 
(repelled  from  the  crystal). 
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In  devices  employing  an  acceleration  grid,  the 


primary  current  Jp  is  always  less  than  the  MCP  output  current  JQ  because 


the  grid  intercepts  a  portion  of  the  MCP  output  current.  If  n( I  ,V.  ,V  ,V  ) 

W  K  d  y 

is  the  primary  electron  energy  distribution  at  the  surface  of  the  crystal, 
then 


V 


'm 


e"<Iw’VVVdVk 


(4) 


Additionally,  an  effective  collected  secondary  electron  emission  ratio  J'  may 
also  be  defined  by 


A, 


i 


5'<VVV 


e5'fVk-Va*Vg)n<Iw*Vk*Va-Vg)dVk 


(5) 


where  5’*(v|c»va»Vg)  is  the  ratio  of  secondary  electrons  collected  by  parts 
of  the  system  other  than  the  crystal  to  monoenergetic  primary  electrons 
with  MCP  ejection  energy  eV^.  Note  that  S’"  will  be  greater  or  smaller 
than  unity  depending  through  Eq . (5 )  on  Vk ,Va ,Vg ) ,  the  primary  electron 
distribution  n(Iw>v|c>Va,Vg) >  the  MCP  bias  voltage  -V^.and  the  crystal  bias 
voltage  which  is  used  to  control  V^.  Further  details  of  the  secondary 
electron  emission  process  can  be  found  in  Ref. 2  (Appendix  I). 

It  is  the  gap  current  density  Jg  which  determines  the  rates  at  which 
electrons  accumulate  on  or  are  removed  from  the  dielectric  mirror  surface. 


Jg  is  given  by 


jg  =  Jpo-n. 


(6) 


Figure  3  is  a  sketch  of  the  gap  current  density  as  a  function  of 


Vg  ,  with  Iw>Va  and  Vm  constant,  for  a  device  employing  a  grid  and  operated 


such  that  V  „  is  positive, 
eq 


v  >  :  a  :  : : 


Notice  that  there  are  four  regions  of  operation  in  Fig.  3.  In  region  4, 
for  large  positive  values  of  Vg  (achieved  with  large  V^),  7'  is  always 
less  than  unity  because  most  of  the  secondary  electrons  generated  are  attracted 
back  to  the  dielectric  mirror.  In  region  3  (the  range  spanning  small  positive 
and  small  negative  values  of  Vg)  most  of  the  secondary  electrons  generated 
will  escape  from  the  dielectric  mirror  surface  and  7"  can  be  greater  than 
unity.  For  more  negative  values  of  Vg  ,  there  is  a  region  (region  2  )  in 
which  the  energy  of  the  majority  of  the  primary  electrons  is  so  low  that 
very  few  secondary  electrons  are  generated  and  thus  7'  is  again  less  than 
unity.  Finally,  in  region  1  ((Vg  <  -VL)  where  VL  =  Va  +  Vk(max),  and  eVk  (max) 
is  the  ejection  energy  of  the  most  energetic  MCP  output  electron)  all  the  primary 
electrons  are  repelled  from  the  dielectric  mirror  surface  and  the  gap  current 
is  then  identically  zero.  This  condition  is  referred  to  as  lockout. 

The  point  E  (corresponding  to  7"  =  1)  at  the  boundary  between  regions 
4  and  3  in  Fig.  3  is  a  stable  operating  point  because  (a)  in  region  4,  as 
electrons  accumulate  on  the  mirror,  Vg  falls,  thereby  decreasing  the  attrac¬ 
tion  of  secondary  electrons  to  the  mirror  and  thus  7'  increases  towards 
unity,  and  (b)  in  region  3  as  the  mirror  loses  electrons,  Vg  increases  and 
attracts  a  larger  fraction  of  the  secondary  electrons  back  to  the  surface 
of  the  mirror,  and  so  7'  decreases  towards  unity.  By  similar  arguments,  it 
can  be  seen  that  the  point  F  in  Fig.  3  is  an  unstable  equilibrium  point. 

During  the  standard  write  and  erase  processes,  the  lockout  regime  is  avoid¬ 
ed  and  the  device  is  operated  in  regions  3  and  4.  However. in  the  space-domain 
or  internal  processing  mode,  sophisticated  operations  such  as  multilevel  thresh¬ 
olding  is  achieved  by  driving  certain  parts  of  the  image  into  the  lockout 
regime  during  operation. 
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WRITE  MODES 

As  illustrated  in  Fig.  3,  may  be  either  positive  or  negative 
depending  on  the  operating  conditions  and,  therefore,  positive  or  negative 
surface  charge  distributions  may  be  written  on  the  dielectric  mirror. 

Thus  two  write  modes  are  possible:  the  electron  accumulation  mode  in 
which  the  mirror  surface  gains  the  appropriate  electron  distribution,  and 
the  electron  depletion  (or  secondary  electron  emission)  write  mode  in  which 
an  electron  distribution  is  extracted  from  the  surface  of  the  mirror. 


Electron  Accumulation  Write  Mode 


In  the  electron  accumulation  write  mode,  which  is  illustrated 

in  Fig.  4,  Vfa  initially  divides  capacitively  between  the  crystal  and 

the  gap  [Fig.  4(a)]  thereby  raising  Vg  into  region  4  of  Fig.  3  for  which 

Z'  <  1.  Then  when  the  write  light  is  turned  on,  electrons  are  deposited  on 

the  crystal,  V  begins  to  fall  toward  its  equilibrium  value  V  ,  and  the 
9 

crystal  voltage  Vx  rises  as  shown  in  Fig.  4(b).  Electron  accumulation  may 
be  terminated  when  the  desired  Vx  is  reached  by  turning  off  the  write  light. 
To  avoid  saturation,  no  parts  of  the  image  must  be  allowed  to  reach  equili¬ 
brium.  The  crystal  may  be  erased  via  secondary  electron  emission  by  flood¬ 
ing  the  photocathode  with  light  and  lowering  Vg  (by  lowering  Vb)  into 
region  3  of  Fig._3,  where  Z'  >  1  [see  Fig.  4(c)].  In  practice,  to  avoid 
lockout,  Vb  is  ramped  downward  no  faster  than  is  driven  to  Vgq  by  the 
electron  removal  process. 
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Electron  Depletion  Write  Mode 


In  the  electron  depletion  write  mode,  which  is  illustrated  in  Fig.  5, 
the  grid  voltage  Va  is  set  at  some  positive  value  which  may  be  several  kilo- 

a 

volts  depending  on  the  halfwave  voltage  of  the  crystal.  Vb  is  then  set  at 
some  initial  bias  level  Vbb  and  the  photocathode  illuminated  uniformly  with 
light  so  that  electrons  are  deposited  or  removed  from  the  mirror  until  all 
points  on  the  mirror  surface  are  driven  to  equilibrium  and  the  desired 
readout  background  level  [see  Fig.  5(a)]  is  obtained. 

With  the  write  information  incident  on  the  photocathode,  Vfa  is  abruptly 
switched  from  the  bias  value  Vbb  to  the  write  value  Vbw,  so  that  Vg  for  all 
points  on  the  crystal  shift  into  region  3  of  Fig. 3  where  7'>1.  Electron  remov¬ 
al  then  begins  at  a  rate  which  is  initially  proportional  to  the  write  light 
intensity  Iw  [see  Fig.  5(b)]. 

Unfortunately,  since  0  is  a  function  of  V  ,  the  electron  removal  rate 

y  y 

will  not  continue  to  be  exactly  proportional  to  I  throughout  the  write 

n 

process.  A  more  sophisticated  control  of  Vb,  however,  can  minimize  this 
problem.6  To  erase  the  image,  Vx  can  be  restored  to  its  original  uniform 
level  by  again  raising  Vb  up  to  Vbb  and  flooding  the  photocathode  with  light 
[see  Fig.  5(c)].  In  practice,  the  value  of  Vg  is  chosen  such  that  AVb  (where 
AVb  “  ^bb  "  ^bw^  ^ar9er  than  or  e£)ua^  t0  halfwave  voltage  of  the  crystal 
if  full  modulation  depth  is  desired.  Because  the  device  is  operated  so  that  no 
areas  of  the  crystal  become  locked  out,  this  form  of  operation  is  called  the 
non-lockout  mode. 


Fig  5.  Illustrating  the  dynamic  characteristics  of  the  electron  depletion 
write  mode,  (a)  Vb  set  at  Vhb  and  mirror  surface  written  to  equili 
brium  to  give  a  uniform  background;  (b)  V5  switched  to  Vbw»  and  Vg 
approaching  Veq  with  the  wright  light  on;  (c)  Vb  switched  back  tos 
Vbb  and  the  image  erased  by  electron  accumulation. 
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READOUT  MODES 

Optical  information  that  has  been  written  into  the  MSLM  may  be  read  out 
with  crossed  polarizers,  interferometrically,  or  with  Schlieren  optics.  The 
possible  readout  modes  depend  on  the  symmetry  and  cut  of  the  electrooptic  crystal 
For  amplitude-only  modulation,  polarization  readout  may  be  used  in 
conjunction  with  a  crystal  cut  for  which  the  components  of  the  allowed 
polarizations  in  the  crystal  experience  phase  changes  with  equal  magnitudes 
but  opposite  sign.  For  phase-only  modulation,  polarization  readout  may  be 
employed  with  the  polarizer  and  analyzer  oriented  along  the  same 
induced  dielectric  axis  of  the  crystal.  However,  for  materials  in  which 
the  allowed  polarization  components  experience  identical  phase  changes 
(e.g.  in  z-cut  3m  crystals  with  longitudinal  fields  only)  it  would  not  be 
necessary  to  polarize  the  readout  light.  Interferometric  and  Schlieren 
readout  schemes  always  lead  to  simultaneous  amplitude  and  phase  modulation 
(intensity  modulation). 

The  readout  mode  that  was  used  in  the  early  prototype  devices  employed 
the  interference  between  the  front  and  back  surface  reflections  from  a 
plane-parallel  electro-optic  crystal.2  If  dielectric  mirrors  are  deposited 
on  both  sides  of  such  a  plane-parallel  crystal  to  form  a  Fabry-Perot 
etalon  of  finesse  F,  a  complete  off-to-on  transition  can  be  achieved  with 
approximately  1/F  times  the  halfwave  surface  charge  density  of  the  crystal, 
and  a  concommitant  increase  in  framing  speed  results.  The  sensitometry  of 
such  a  device  is  distinguished  by  a  very  large  value  of  ganrnia,  where  y  is 
defined  by 

y  *  iog(i01/i02)/iog(E2/E1)  (7) 

Here,  E^  and  Eg  and  Ig^  and  Igg  are  the  corresponding  exposures  and  read¬ 
out  intensities  for  two  points  on  the  linear  portion  of  the  device 
sensitometry  curve.  Figure  6  is  a  plot  of  the  readout  intensity  vs.  surface 


charge  density  of  an  MSLM  with  a  Fabry-Perot  crystal  etalon  for  two  different 
surface  reflectivities. 

CYCLING  MODES  AND  OTHER  ERASE  MECHANISMS 

The  MSLM  can  be  either  framed  through  a  discrete  sequence  of  write, 
read,  and  erase  periods  or  operated  in  a  continuous  mode.  In  the  framed 
mode,  the  optical  write  signal  is  integrated  and  stored  in  the  form  of  a 
charge  distribution  on  a  high  resistivity  surface;  after  a  readout  period  of 
the  desired  length  has  elapsed,  the  charge  distribution  is  removed- 

In  addition  to  the  erase  methods  based  on  electron  accumulation  and 
depletion  using  the  MCP  and  photocathode  as  the  source  of  electrons,  several 
other  erase  mechanisms  can  be  employed.  For  example,  high  resistivity 
crystals  with  high  resistivity  dielectric  mirrors  may  be  coated  with  an 
upper  layer  of  photoconductive  material  which  has  a  high  dark  resistivity. 

If  a  ring  electrode  is  placed  on  the  photoconductor  near  the  edge  of  thisstruc 
ture,  the  charge  distribution  can  be  removed  by  flooding  the  photocon¬ 
ductive  layer  with  light  and  appropriately  biasing  the  ring  electrode. 

The  continuous  mode  of  operation  can  be 

achieved  with  low  resistivity  and  photoconductive  crystals  by  allowing  the 
charges  to  bleed  through  the  crystal.  In  the  continuous  mode,  the  instan¬ 
taneous  modulation  is  proportional  to  the  write  image  intensity  for  all 
temporal  variations  within  the  bandwidth  of  the  system.  However,  large  MCP 
output  currents  are  required  to  operate  the  device  in  the  continuous  mode,  so 
this  mode  of  operation  is  not  recommended  for  high  framing  speed  applications. 


SPACE-DOMAIN  PROCESSING  MODE 

By  exploiting  the  secondary-electron-emission  characteristics  of  the 
dielectric  mirror  and  the  charge-readout  intensity  characteristic  of  the 
crystal,  the  MSLM  can  be  operated  in  its  space-domain  information  processing 
mode.  In  this  internal  processing  mode,  the  device  can  perform  several 
sophisticated  information  processing  operations  on  information  that  has 
already  been  stored  in  the  device  and  on  information  as  it  is  being  read 
(in  real  time)  into  the  device. 

The  space-domain  nonlinear  operations  of  analog  thresholding,  real¬ 
time  hardclipping,  and  contrast  reversal  are  the  key  to  several  of  the 
unique  space-domain  information  processing  operations  that  are  achievable 
with  the  MSLM.  These  innate  operations  of  the  MSLM  are  a  consequence  of 
the  secondary  electron  emission  processes  that  are  associated  with  the 
electron  beam  addressing  of  the  floating  surface  of  the  dielectric  mirror. 

Figure  7  is  a  graphic  illustration  of  the  above-mentioned  nonlinear 
operations  as  they  are  achieved  by  the  MSLM.  By  combining  these  three  opera¬ 
tions  with  the  standard  spatial  light  modulator  operations  of  (1)  multiplication 
of  the  write  and  readout  beam  intensities,  and  (2)  addition  of  two  or  more 
write  beams,  plus  the  fact  that  multiple  it  radians  of  phase  retardation  are 
achievable  with  proper  choice  of  crystal,  the  MSLM  has  been  able  to  perform 
space-domain  image  subtraction,  stored-image  hard  clipping,  real-time  and 
stored-image  contrast  enhancement,  real-time  and  stored-image  edge  enhancement, 
as  well  as  binary-level  optical  logic  operations  such  as  AND,  OR,  NAND,  NOR, 

XOR  and  NXOR.  The  principles  and  the  limitations  of  the  three  key  operations 


are  described  below. 


Graphical  Illustration  of  the  cross  sectional  profile  of  (a)  ar,  ori¬ 
ginal  Image;  (b)  the  Image  after  the  MSLM  analog  thresholding  opera¬ 
tion;  (c)  the  Image  after  hard  clipping;  (d)  the  image  after  contras 
reversal.  Ij  Is  the  threshold  Intensity. 


Analog  Thresholding 

Consider  two  areas  A'  and  A"  within  an  electron  distribution  on  the 
surface  of  the  crystal  with  charge  densities  o'  and  c" ,  respectively.  If 
o"  is  sufficiently  greater  than  o',  then  can  be  adjusted  so  that,  when 
the  photocathode  is  uniformly  illuminated,  primary  electrons  from  the  MCP 
are  repelled  from  A'  but  not  from  A".  Under  these  conditions  A"  and  all 
areas  with  charge  densities  greater  than  a"  will  be  erased,  while  A'  and  all 
other  areas  for  which  o<a‘  will  be  unaffected.  This  operation  (analog 
level  thresholding)  corresponds  to  erasing  all  the  parts  of  an  image  for 
which  the  intensity  falls  below  some  threshold  intensity  Iy  as  illustrated 
graphically  in  Fig.  7(b).  That  is. 


Iw(x,y)  — *  Io(x'y)  =  A^w,max  ’  Iw(x’y)]  <8> 

where  I  is  the  maximum  value  of  Iw(x,y),  the  write  light  intensity, 

and  A  is  a  constant.  Repeated  applications  of  this  operation  (multiple-level 
thresholding)  can  be  used  to  analyze  the  spatial  distribution  of  intensities 
in  an  image.  Analog  thresholding  is  useful  in  the  implementation  of  contrast 
enhancement  and  the  logic  operations  AND  and  NAND. 

Real-Time  Hard  Clipping 

This  operation  is  most  conveniently  performed  in  the  electron  depletion 
(secondary  emission)  write  mode.  The  MSLM  is  initially  uniformly  biased  in 
the  OFF  state  with  charge  density  Og.  Then  Vb  is  ramped  downward  (from 
higher  to  lower  voltage)  at  a  rate  with  the  optical  image  incident  on 
the  photocathode.  All  the  intensities  below  some  threshold  level  IT  (which 
depends  on  Vb)  will  be  barely  recorded  because,  at  their  locations,  electrons 
cannot  be  removed  fast  enough  from  the  crystal  to  prevent  the  gap  field 
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from  eventually  repelling  all  future  primary  electrons.  Thus,  all  intensi¬ 
ties  below  threshold  will  remain  at  a  low  intensity  level  (S  =  0  state). 

Because  the  slope  of  the  Vx  vs  write  exposure  characteristic  can  be  very  steep 
between  threshold  and  saturation,  virtually  all  write-light  intensities 
above  threshold  are  written  to  the  maximum  charge  density  which  is  set  to 
coincide  with  the  S  =  1  state.  The  hard-clipping  operation  is  illustrated 
graphically  in  Fig.  7(c),  where 


Iw(x,y) 


I0(x,y)  = 


'1  Iw(x,y)  >  IT 

0  IJx.y)  <  I. 


(9) 


‘w'"’*7'  Aj 

This  operation  can  lead  to  real-time  contrast  enhancement  if  the  crys¬ 
tal  has  a  low  halfwave  voltage  (e.g.  DKDP  at  -51°C,  VTr<  200V)  and  if  Vb 
is  adjusted  so  that  the  spatial  intensity  variations  of  the  image  are  mapped 

onto  the  knee  of  the  Vv  vs  Iw  curve.6 

X  w 


Contrast  Reversal 

There  are  two  methods  of  performing  contrast  reversal.  Since  the 
two  write  modes  permit  images  to  be  written  as  negative  or  positive  charge 
distributions,  pockels-effect-type  crystals,  because  they  exhibit  the 
linear  electro-optic  effect,  will  yield  contrast  reversed  images  when 
operated  in  the  accumulation  and  electron  depletion  write  modes.  If  the 
image  is  already  stored  in  the  MSLM,  then  because  the  readout  intensity 
characteristic  is  periodic  in  surface  charge  density  under  polarization 
and  interferometric  readout,  contrast  reversal  can  be  achieved  by  adding 
the  appropriate  uniform  charge  density  to  dielectric  mirror.  In  this 
operation, 

Iw(x,y)  — *  I0(x,y)  =  A[Iw(max)  -  Iw(x,y)],  (10) 

where  I  (max)  is  the  maximum  value  of  I  (x,y),  the  write  light  intensity, 
and  A  is  a  constant.  This  operation  is  illustrated  graphically  in  Fig.  7(d). 


Figure  8  suiranarizes  some  of  the  results  of  the  space-domain 
information  processing  operations  achievable  with  the  MSLM. 

Figure  8(a)  shows  the  image  of  an  U.S.  Air  Force  resolution  chart. 

The  sequence  of  Figs .8(b)-8(d)  illustrates  analog  and  digital  space-domain 
intensity  thresholding.  Figure  8(b)  shows  a  four-level  image.  Figure  8(c) 
shows  the  same  image  after  space-domain  thresholding  in  which  quadrant  1 
was  reduced  to  the  zero-intensity  level  of  quadrant  3,  leaving  quadrants  2 
and  4  unchanged.  Figure  8(d)  is  the  resultant  binary-level  image  after  the 
MSLM  with  a  stored  image  of  Fig. 8(c)  was  illuminated  uniformly  with  light 
to  bring  quadrant  4  to  the  level  of  quadrant  2  and  thresholding  was  again 
employed  to  restore  quadrants  1  and  3  back  to  zero. 

Figures  8(e)  and  (f)  illustrate  contrast  reversal.  In  both  cases,  the 
images  were  written  by  electron  accumulation.  The  crystal  was  first  biased, 
by  means  of  secondary  electron  emission,  with  a  uniform  positive-charge 
background  at  either  a  minimum  in  the  transmittance  characteristic  [Fig. 8(e)] 
or  at  maximum  transmittance  [Fig. 8(f)]. 

Figures  8(g)  and  8(h)  illustrate  contrast  enhancement.  In  Fig. 8(g) 
the  image  shown  in  Fig. 8(e)  was  washed  out  with  background  illumination. 
Figure  8(h)  shows  the  contrast-enhanced  result  after  analog  thresholding. 
Figures  8(i)  and  8(j)  illustrate  edge  enhancement  of  a  binary-level  image. 
Figure  8(i)  results  from  Fig. 8(e)  when  the  write-image  exposure  corresponds 
to  the  full-wave  surface-charge  density  02^;  similarly.  Fig. 8 ( j )  results 
from  Fig. 8(f)  for  the  full-wave  input-image  exposure. 

Figures  8(m)-8(r)  illustrate  the  AND,  NAND,  OR,  NOR,  XOR,  and  NXOR 
operations  of  the  horizontal  and  vertical  bar  images  shown  in  Figs.8(k)  and 
8(it).  Except  for  Figs.8(o)  and  8(p).  all  the  images  were  achieved  in  a 
device  without  an  acceleration  grid.  Further  details  of  how  the  above 


Fig  8.  Some  of  the  open-loop  information  processing  operations  that  have 
been  achieved  by  operating  the  MSLM  in  its  space  domain  informati 
processing  mode. 


operations  were  achieved  are  given  in  Ref.  2. 

The  image-processing  operations  presented  in  Fig. 8  were  obtained  using 
a  vacuum-demountable  prototype  device  with  a  330-um-thick  oblique-cut  Li Nb03 
crystal.  This  device  employed  a  single  dielectric  mirror  on  the  crystal, 
and  was  read  out  between  crossed  polarizers.  The  MCP  used  had  a  strip 
current  of  250yA  and  a  gain  of  1.1X104  at  IkV  of  bias.  The  device  was 
studied  with  and  without  an  acceleration  grid.  The  grid  had  a  spatial 
frequency  of  30  wires /mm. 

The  strip-current-limited  framing  rate  of  this  device  with  the  grid  at 
2  kV  was  30Hz  at  full  modulation  depth  and  60Hz  at  50%  modulation  depth.  The 
corresponding  halfwave  write  (electron  accumulation)  and  erase  times  were 
10  and  23  msec,  respectively.  Its  measured  normalized  light-to-light  coher¬ 
ent  MTF  was  found  to  have  a  value  of  approximately  1.9  cycles/mm  at  50%  con¬ 
trast  for  a  gap  £g  of  500ym.  Similar  resolution  results  were  obtained  with 
the  gridless  device. 

Limitations  of  the  Space-domain  Processing  Mode 

The  fidelity  of  the  space  domain  processing  operations  described  above  is 
strongly  dependent  on  the  energy  distribution  n(Iw,  Vk,  Vfl,  V  )  of  the  primary 
electrons,  the  secondary  electron  emission  characteristics  of  the  grid  and  the 
surface  that  receives  the  array  of  primary  electron  beams,  the  readout-intensity- 

vs. -surface-charge-density  characteristic,  and  the  spatial  resolution  of  the 
device. 


Sharp  well-defined  thresholds  are  achieved  when  the  threshold  volt¬ 
age  spread  AVy  is  much  smaller  than  the  crystal  halfwave  voltage  V,R. 
Ideally,  for  a  sharp,  well-defined  thresholding  level,  the  array  of 
primary  electron  beams  incident  on  the  mirror  should  be  monoenergetic. 
Additionally,  the  secondary  electron  emission  ratio  for  the  grid  should 
be  small  and  that  for  the  mirror  surface  should  have  a  maximum  value  that 
is  much  greater  than  unity.  Unfortunately,  the  majority  of  the  output 
electrons  from  straight-channel  microchannel  plates  have  a  considerable 
energy  spread  (Oto > 50eV) ,s  Thus,  sharp  thresholds  are  not  expected  from 
straight-channel  MCPs.  Better  defined  threshold  levels  are  expected  from 
curved  channel  MCPs  which  should  have  a  much  smaller  output  electron 
energy  spread.  In  the  current  oblique-cut  LiNb03  MSLMs  that  employ  straight 
channel  MCPs,  AVT  s  lOOV  and  V  =  1250V. 

The  readout-intensity-vs-surface-charge-density  characteristic  can 
also  be  of  considerable  importance  in  space  domain  processing  when  the 
electron  accumulation  write  mode  is  used.  For  example,  the  sharp  OFF  to 
ON  characteristic  of  the  Fabry-Perot  crystal  [see  Fig.  (6)]  permits  hard 
clipping  to  be  accomplished  in  the  electron  accumulation  write  mode  when 
the  initial  readout  background  or  bias  intensity  is  set  at  a  null. 

Severe  degradation  in  image  sharpness  can  result  during  space  domain 
processing  if  care  is  not  exercised  to  prevent  stray  secondary  electrons 
from  being  attracted  back  to  the  mirror  surface,  or  if  the  gap  voltage  Vg 
becomes  so  small  that  proximity  focussing  of  the  primary  electrons  is  lost. 
It  is,  therefore,  advantageous  to  have  a  priori  knowledge  of  the  intensity 
dynamic  range  of  the  write  information  when  real-time  space  domain  proces¬ 
sing  is  being  attempted. 
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MATERIALS  LIMITATIONS 

The  versatility  of  the  MSLM  is  derived  from  the  unique  architecture 
of  the  device  which  permits  one  to  choose  materials  for  the  photocathode,  MCP, 
grid,  dielectric  mirror  and  the  electro-optic  crystal  that  are  best  suited  for 
the  application  in  question.  We  shall  discuss  each  of  these  in  turn. 
Photocathode 

The  input  window  and  photocathode  material  determine  the  write  light  wave¬ 
length  response  of  the  MSLM.  Typical  photocathode/window  combinations 

span  the  wavelength  range  between  110  and  1200  nm.  Fortunately,  MCPs  respond 
directly  to  radiation  in  the  mid  ultraviolet  and  soft  x-ray  regions  and  there¬ 
fore  no  photocathode  is  needed  in  these  regions.7’ 3  Beyond  1300nm  no  photo¬ 
cathode  exists  and  consequently  the  MSLM  cannot  process  infrared  signals  with 
wavelengths  beyond  this  cutoff  value. 

MicroChannel  Plate 

A  microchannel  plate  (MCP)  is  an  imaging  electron  multiplier  that  preserves 
the  relative  spatial  distribution  of  the  input  electron  flux.  It  is  essentially 
a  thin  (0.5  -  1.0  mm)  wafer  of  semiconducting  glass  channels  that  is  electroded 
on  both  faces.  A  typical  MCP  may  have  10-um-diameter  pores  on  15  ym  center  to 
center  spacing  packed  in  a  hexagonal  array.  Both  straight  and  curved  channel9 
MCPs  are  manufactured.  In  straight  channel  plates,  the  pores  are  generally 
slanted  at  an  angle  of  about  5-15°  away  from  the  plate  normal.  The  typical 
length-to-diameter  ratio  of  the  pores  is  about  40  to  1.  Plates  as  large  as 
15cm  in  diameter  have  been  manufactured. 

The  strip  current  I$  is  that  current  which  flows  through  the  plate  when 
no  electron  or  photon  is  incident  on  the  plate.  The  typical  strip  current  for  a 
25-mm-diameter  straight  channel  plate  is  about  10yA  at  IkVDC  of  bias.  Under 
steady  state  dc  operation,  the  unsaturated  MCP  output  current  density  J  is 


Jo  =  !w  neG/hv 


where  I  and  v  are  the  intensity  and  frequency  of  the  write  light  respec- 

W 

tively,  h  is  Planck's  constant,  e  is  the  electronic  charge,  n  is  the  quantum 
efficiency  of  the  photocathode  and  G  is  the  gain  of  the  MCP.  In  continuous 
mode  operation  the  extracted  current  from  the  MCP  is  generally  less  than  one- 
fifth  of  the  strip  current. 

In  pulsed-mode  applications,  pulsing  the  strip  current  can  result  in  out- 
put  current  pulses  as  high  as  1  A/cm  . 10  For  high  gain  applications,  two 
straight  channel  plates  may  be  cascaded  in  a  proximity  focus  configuration. 
Recently,  curved  channel  plates3  operating  at  2  kV  of  bias  with  strip  currents 
of  30  uA  and  gains  of  106  have  been  replacing  the  cascaded  configuration. 

The  framing  speed  of  the  MSLM  is  generally  limited  by  the  available  out¬ 
put  current  of  the  MCP.  And,  since  the  MCP  output  current  is  derived  from  the 
strip  current,  high  strip  current  is  often  more  desirable  than  high  gain.  The 
MCP  strip  current  depends  on  the  resistivity  of  the  channels,  and  special  post¬ 
manufacturing  reduction  processes  can  be  employed  to  lower  the  resistivity  of 
the  pores  and  thereby  increase  the  strip  current  without  a  substantial  reduc¬ 
tion  of  the  gain  of  the  plate.  Straight  channel  MCPs  with  strip  currents  up  to 
300  yA  have  been  employed  in  MSLMs.  To  avoid  thermal  runaway,  care  must  be 
exercised  when  using  plates  with  such  high  strip  currents. 

Since  there  has  not  been  a  great  demand  for  high  strip  current  MCP's,  very 
few  of  them  have  been  produced,  and  little,  if  any,  systematic  research  has  been 
done  on  the  development  of  these  items.  Such  an  effort  is  clearly  warranted  at 


Acceleration  Grid 


The  ideal  requirements  of  the  acceleration  grid  are  large  fractional 
open  area,  high  tensile  strength,  low  elasticity,  and  low  secondary  elec¬ 
tron  emission  yield.  In  general,  a  metal  mesh  (e.g.  of  copper  or  nickel) 
is  chosen  as  the  grid  material. 

High  tensile  strength  is  required  because  the  electrostatic  forces 
acting  on  the  grid  during  operation  can  be  quite  large.  More  specifically, 
it  can  be  shown  that  the  electrostatic  force  per  unit  area  T  acting  on  a 
planar  grid  which  is  held  at  a  potential  Va  and  positioned  at  distances  i  and 
jig  from  the  MCP  output  electrode  and  the  free  surface  of  the  crystal  res¬ 
pectively,  is  given  by11 


F  -  H0[<VV2  -  'W2] 


f^b  -  V  +  °'*0 

2 

' 

Va 

2 

1  + 

.  , 

*m 

k  4 

where  &x  and  eB  are  the  thickness  and  dielectric  constants  respectively  of 
the  crystal  in  the  longitudinal  direction,  and  eg  is  the  permittivity  of 
free  space.  In  Eq.  12,  T  is  positive  when  the  force  is  in  the  direction 
of  the  crystal  and  it  increases  when  electrons  are  removed  from  the  mirror. 

For  typical  MSLM  parameters  Um  -  i  =  0.5mm,  z%  =  0.1mm,  a  =  =  +0.42uC/cm2, 

e()  *  38,  =  5kV,  and  =  2kV)  we  find  F  =  0.025  Newton/cm^. 

Low  elasticity  and  high  tensile  strength  are  therefore  important  require¬ 


ments  so  that  the  grid  does  not  deform  significantly,  which  wqj Id  cause  an  increase  in 
the  force  Clearly,  the  need  to  maintain  a  small  gap  between  the  MCP  and  the 


crystal  to  assure  good  proximity  focussing,  and  at  the  same  time  the  need 
to  employ  a  large  enough  gap  so  that  the  grid  does  not  touch  or  arc  to  either 
surface  of  the  MCP  or  the  crystal  during  operation  are  the  constraints  that 
set  the  upper  and  lower  limit  on  the  dimensions  of  the  gap. 

Dielectric  Mirror 

Dielectric  mirrors  consist  of  alternating  layers  of  high  and  low  re¬ 
fractive  index  dielectric  materials.  In  the  MSLM  design  shown  in  Fig.l, 
the  electrons  are  deposited  on  the  dielectric  mirror.  Ideally,  then,  the 
dielectric  mirror  should  have  high  damage  resistance  under  electron  beam 
bombardment,  negligible  electron  beam  induced  conductivity,  large  secondary 
electron  yield,  high  reflectivity  at  the  readout  wavelength,  low  vapor 
pressure,  and,  to  avoid  cracking  and  flaking,  low  internal  stresses  when 
cycled  up  and  down  in  temperature.  Several  combinations  of  materials  exist 
that  satisfy  these  requirements. 


The  performance  of  the  MSLM  is  determined  more  by  the  electro-optic 
crystal  than  by  any  of  the  other  components.  The  crystal  in  the  MSLM,  like 
that  in  other  electro-optic  spatial  light  modulators,  is  operated  in  the 
longitudinal  electro-optic  mode.  That  is,  the  readout  light  and  the  applied 
electric  field  component  of  interest  are  parallel.  The  electro-optic 
crystal  within  the  MSLM  is  addressed  through  the  dielectric  mirror  by  the 
charge  distribution  on  the  surface  of  the  mirror,  and  therefore,  the  half¬ 
wave  surface  charge  density  is  often  a  better  figure  of  merit  for  comparing 
crystals  than  the  halfwave  voltage. 

The  halfwave  longitudinal  surface  charge  density  is  the  charge 
density  at  the  surface  of  the  dielectric  mirror  that  gives  rise  to  a  phase 
retardation  of  ir  radians  in  the  electro-optic  crystal.  The  halfwave  surface 
charge  density  is  related  to  the  halfwave  voltage  via  the  capacitance 
per  unit  area  C  of  the  crystal  by  a^CV^. 

When  amplitude  and  phase  modulation  are  achieved  by  modulating  the 
birefringence  of  the  crystal  and  reading  out  the  crystal  between  crossed 
polarizers,  the  phase  retardation  for  amplitude  modulation  ra  is  given  by 


r  =  A$v-  -  A <t>v^,  where  A4>v.  and  A<tu-  are  the  phase  changes  produced  along  the 
induced  dielectric  axes,  x'  and  y',  of  the  crystal.  Phase-only  modulation 
with  such  a  crystal  is  achieved  by  reading  out  the  device  with  light  polar¬ 
ized  along  either  the  x'  or  y'  axis.  Then  the  phase  retardation  for  phase 
modulation  rp  is  given  either  by  rp  =  A^'  or  rP  a  In  z~cut  ^2m 

crystals,  for  example,  A$x  *  -A<£y,  hence  Tp  =  ira,  and  so  =  20^.  Also, 
since  the  MSLM  is  a  reflex  mode  device,  the  readout  light  makes  a  double¬ 
pass  through  the  crystal  and  hence  =  icT,  where  o^  is  the  reflex  half¬ 
wave  surface  charge  density. 

Lithium  Niobate 

LiNb03,  a  crystal  of  the  3m  point  group,  has  been  the  electro¬ 
optic  material  most  frequently  used  to  date  in  the  MSLM.  The  symmetry  of 
this  crystal  makes  the  z-cut  material  fairly  well  suited  for  phase  modula¬ 
tion  because  the  longitudinal  component  of  the  electric  field  modulates  A£x. 
and  A>?y>  symmetrically.  That  is,  the  index  ellipsoid  remains  cylindri- 
cally  symmetric  when  longitudinal  fields  are  applied  along  the  z-axis. 
Therefore,  in  principle,  there  is  no  need  to  polarize  the  readout  light 
for  phase  modulation  with  this  crystal  cut.  In  practice,  however,  trans¬ 
verse  fields,  which  are  always  present  in  non-uniform  surface  charge 
distributions,  cause  an  additive  phase  modulation  error.  Fortunately,  the 

dielectric  mirror  has  a  shielding  effect,  since  the  strongest  transverse 
fields  are  confined  within  the  dielectric  mirror  layer. 

To  achieve  intensity  modulation  with  the  z-cut  crystal,  an  interfero¬ 
metric  or  a  Schlieren  readout  configuration  may  be  used.  An  interferometric 


readout  scheme  that  works  fairly  well  involves  interference  between  the 
back  and  front  surface  reflections  (15%  reflectivity  at  normal  incidence) 
from  a  plane  parallel  crystal. 


For  z-cut  LiNb03  operated  in  the  reflex  mode  with  longitudinal  fields, 
the  phase  retardation  rpR  due  to  a  longitudinal  field  Ez  is  given  by 

rpR  ’  4*x  '  %  -  (13) 


and  hence 


vttR  =  2Fr— 

‘•"o'  13 


°ttR  = 


eoeitA 


•o'  13*X 


(14a)  ;  (14b) 


where  e0  is  the  permittivity  of  free  space,  e(|  is  the  dielectric  constant 
of  LiNbOg  along  the  field  direction,  \  is  the  wavelength  of  the  readout 
light,  and  n0,£x  and  r  are  respectively  the  ordinary  refractive  index, 
thickness  and  electro-optic  coefficient  of  the  LiNbC^  wafer.  For  z-cut  LiNbOj, 
rl3=8.6X10"12m/V,  na=2.28  at  633nm  wavelength,  and  £„  =30.  Thus  V7rp=3100  volts, 
and  for  a  lOOynrthick  crystal,  o^O.SyC/cm2 . 

The  halfwave  voltage  and  the  halfwave  surface  charge  density  for  inten¬ 
sity  modulation  can  be  reduced  significantly  by  employing  the  crystal  cut 
which  optimizes  the  phase  retardation  .  It  can  be  shown  that  the  optimum 
cut  is  such  that  the  normal  to  the  plate  lies  in  the  x-y  plane,  and  is 
rotated  about  the  x-axis  by  55°  from  the  +z  towards  the  -y  axis.12  For  this 
cut,  the  caluclated  value  of  is  1150  volts  which  compares  favorably  with 
our  measured  values  of  1250  volts.3  It  should  be  pointed  out  that  for  this 


cut  | A<l>x  J ^ | A4>y  J  ,  and  therefore,  even  with  longitudinal  fields  only,  this 
crystal  cut  will  not  lead  to  an  amplitude-only  modulator.  That  is,  there 
would  always  be  a  signal -dependent  phase  error  accompanying  the  amplitude 
modulation. 


-37- 


Ideally,  for  MSLM  applications  the  electro-optic  crystal  should  have 
small  dielectric  constants,  large  electro-optic  coefficients  which  couple 
only  to  longitudinal  electric  fields,  high  surface  and  bulk  resistivity, 
high  dielectric  strength,  good  optical  quality,  good  mechanical  strength 
and  hardness,  small  temperature  coefficients,  insolubility  in  the  common 
solvents,  low  vapor  pressure  and  good  bakeability  up  to  about  400°C. 

Unfortunately,  because  of  the  coupling  between  the  dielectric  constants 
and  the  electro-optic  coefficients  in  most  materials,  any  steps  taken  to 
increase  the  electro-optic  coefficients  (e.g.  by  operating  near  the  Curie 
temperature)  always  result  in  a  concommitant  increase  in  the  dielectric 
constant,  and,  in  general,  design  tradeoffs  must  be  made.  Additionally, 
both  the  spatial  resolution  and  the  capacitance  increase  with  decreasing 
crystal  thickness,  and  because  the  primary  current  is  limited  by  the  MCP 
strip  current,  tradeoffs  between  speed  and  resolution  are  often  necessary. 
Research  on  new  electro-optic  materials  for  spatial  light  modulation  and  on 
new  ways  of  employing  existing  materials  is  sorely  needed  at  this  time. 


Electro-optic  Materials  Comparison 


Almost  any  electro-optic  material  can  be  employed  in  the  MSLM,  and  the 
choice  of  material  would  be  dictated  by  the  intended  application.  Nevertheless, in 
Table  I  we  list  and  compare  some  of  the  properties  of  importance  to  the 
MSLM  for  the  common  electro-optic  crystals  used  in  spatial  light  modulators 
today.  The  materials  chosen  for  comparison  are  z-cut  deuterated  potassium 
dihydrogen  phosphate  (DKDP),  which  is  used  in  the  photo  DKDP  light  valve,13 
<001>-cut  Bismuth  silicon  oxide  (BSO),  which  is  used  in  PROM,14  and  z-cut 
and  55°-cut  lithium  niobate,  which  have  been  used  in  the  MSLM.2*3  DKDP  is 
generally  used  near  its  Curie  temperature  (-51°C)  in  spatial  light  modulator 
applications. 

Table  I.  Materials  Comparison 
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(25°C) 
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29 

44 

36 
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0.42 

1.2 

0.51 

In  the  table,  e  and  e  ,  and  p  and  p  are  the  dielectric  constants 
ii  -*•  II  * 

and  resistivities  parallel  and  perpendicular  respectively  to  the  longitudi¬ 
nal  electric  field.  To  a  good  approximation,  it  can  be  shown  that  for 
electro-optic  materials  with  large  dielectric  constants  the  MSLM  spatial 
resolution  MTF  rollover  point  (cycles/mm)  is  given  by 


Ni  * i 


Here,  the  quantity  (e^/e^)*  is  used  as  a  crystal  resolution  parameter  in 
Table  I.  Note  that  by  cooling  DKDP  near  its  Curie  temperature,  its 
resolution  parameter  increases  significantly. 

The  halfwave  voltage  and  the  halfwave  surface  charge  density  are  also 
given  in  Table  I.  Note  that  because  z-cut  LiNb03  remains  uniaxial  when  a 
longitudinal  field  is  applied,  we  have  chosen  to  calculate  its  halfwave  voltage 
and  charge  density  that  correspond  to  an  absolute  change  in  phase  of  tt  radians 
for  either  x  or  y  polarized  light.  Note  also  that  55°-cut  LiNb03  offers  the 
lowest  halfwave  surface  charge  density  and  hence  the  reason  for  our  interest 
in  this  crystal  cut.  On  the  other  hand,  the  switching  charge  per  pixel 
which  is  also  provided  in  Table  I,  is  given  by  Q  =c  D/N2  ,  and  is  smallest  for 
for  cooled  DKDP.  The  write  light  switching  energy  per  pixel  7 '  is  given  by 

_  Q  hv 

E,  rar  <16> 

where  we  have  assumed  Jp  =  Jq  .  For  the  55°-cut  LiNbO^  example  shown  in  Table  I, 
taking  A  =  633nm,  n  =  o.l  and  G  =  10^, we  find  that  the  calculated  value  of 
TT  is  l.OpJ/pixel. 


APPLICATIONS  TO  ADAPTIVE  PHASE  COMPENSATION 


Atmospheric  optical  communications  theory15-19  and  propagation  studies20-23 
suggest  that  adaptive  phase  compensation  may  considerably  improve  the  per¬ 
formance  of  optical  communications  links  through  such  low-visibility  conditions 
as  snow,  rain,  fog  or  smoke.  An  adaptive  phase  compensated  receiver  essentially 
senses  the  spatial  phase  distortions  in  the  scattered  light  arriving  at  the 
receiver,  and  uses  a  spatial  phase  modulator  to  "straighten  out"  the  distor¬ 
tions  in  the  wavefront.  The  compensated  wave  can  then  be  processed  by  a 
diffraction-limited  field  of  view  (FOV)  receiver  (e.g.  heterodyne)  with  very 
little  collection  of  the  background  light  which  arrives  from  a  wide  FOV. 
Therefore,  such  a  scheme  results  in  an  improvement  of  the  communication  signal- 
to-noise- ratio. 

A  million  or  more  compensated  elements  are  generally  required  to  correct 
a  low-visibility  receiver  aperture  of  only  a  few  centimeters  diameter.19-23 
The  conventional  adaptive  optical  technology24  consisting  of  an  array  of 
discrete  phase  modulators  cannot  practically  be  extended  to  these  high  resolu¬ 
tions.  Therefore,  an  "all-optical"  parallel  processing  approach,  which  employs 
a  sensitive  optically-addressed  spatial  light  modulator  (MicroChannel  Spatial 
Light  Modulator)  and  a  new  phase  measurement  and  compensation  technique 
(Interference  Phase  Loop),  was  investigated. 

All  Optical  System 

Figure  9  is  a  conceptual  diagram  of  an  "all  optical",  phase-compensated 
direct-detection  communications  receiver.  The  received  wavefront  of  the 
scattered  light  passes  through  an  optically-addressed  monolithic  spatial 
light  modulator  which  adds  the  appropriate  phase  to  present  a  compensated, 
approximately  planar  wavefront  to  the  direct-detection  receiver.  Part  of 
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Fig  9.  Direct-detection  phase  compensated  receiver. 
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Fig  10.  An  all-optical  high  resolution  interference  phase  loop  employing  an 
MSLM. 


the  compensated  wavefront  is  fed  back  into  the  box  labeled  "phase  estimator" 
whose  output  is  an  optical  signal  proportional  to  the  phase  of  the  compensated 
beam.  This  optical  control  signal  drives  the  spatial  light  modulator. 

Most  existing  phase  measurement  techniques5  either  have  undesirable 
performance  limitations  or  are  too  complicated  to  be  implemented  in  the  phase 
estimator  of  Fig.  9.  For  example,  elementary  phase  visualization  techniques 
such  as  Zernike  phase-contrast,  heterodyne,  homodyne,  shearing,  and  polari¬ 
zation  interferometry  are  all-optical;  however,  they  cannot  by  themselves 
distinguish  wavefront  phase  from  amplitude  variations.  They  also  suffer  from 
a  nonlinear  dynamic  range  of  tt  radians  as  well  as  from  2m  and  phase  quadrant 
ambiguities. 

Recently,  we  have  recognized  (see  Refs.  25  and  26  which  are  reproduced 
in  Appendices  III  and  IV)  that  there  exists  a  class  of  closed  loop  configu¬ 
rations  similar  to  Fig.  9  that  possesses  inherent  properties  which  simply 
phase  estimation  and  compensation.  In  these  closed  loop  configurations,  called 
Interference  Phase  Loops  (IPL),  phase  estimation  and  phase  compensation  are 
achieved  simultaneously  by  using  the  output  an  elementary  interferometer  to 
drive  a  phase-only  modulator  in  the  path  of  the  wavefront  to  be  measured  or 
compensated.  Figure  9  therefore  becomes  an  example  of  an  IPL  when  the  "phase 
estimator"  is  an  elementary  two-beam  interferometer. 

Theoretical  and  experimental  studies  have  shown5’25’26  that  the  IPL  not 
only  largely  ignores  amplitude  fluctuations,  but  also  ha*-  ^  phase  quadrant 
ambiguity  and  can  continuously  estimate  phase  over  multiple  tt  radians  of 
dynamic  range.  When  self-interference  is  employed  (e.g.  Zernike  phase  con¬ 
trast),  monochromatic  light  is  nor  required.  In  operation,  the  drive  signal 
to  the  IPL  phase  modulator  and  the  output  intensity  from  the  elementary 
interferometer  become  representations  of  the  phase  estimate.  In  many  respects. 


the  IPL  can  be  considered  as  an  optical  phase-locked  loop  which  directly 
tracks  optical  wavefront  phase. 

Figure  10  illustrates  a  specific  all-optical  implementation  of  the  IPL 
employing  a  homodyne  interferometer  and  a  reflection  readout  spatial  phase 
modulator,  the  MSLM.  Figure  11  is  a  low-resolution  demonstration  of  all- 
optical  phase  compensation  with  an  MSLM/IPL  system.  The  prototype  MSLM  that 
was  used  employed  an  electrooptic  crystal  (LiTa03)  with  parallel  faces  and 
without  a  dielectric  mirror.  The  reflection  from  the  front  crystal  face 
provided  a  homodyne  L.O.  beam,  which  allowed  the  nonuniform  internal  phase 
distortion,  <f>g,(x,y),  of  the  crystal  to  be  compensated.  (This  parallel-face 
structure  is  actually  a  low  finesse  Fabry-Perot  interferometer.6)  Figure  11a 
is  an  i nterferogram  of  the  initially  uncompensated  phase  distribution,  which 
had  a  dynamic  range  of  about  tr  radians.  The  compensated  i nterferogram  in 
Fig.  lib  was  recorded  with  the  same  exposure  as  Fig.  11a,  and  the  longer 
exposure  (12  times  larger)  of  Fig.  11c  reveals  that  the  residual  phase  error 
is  very  small  (^X/10  or  A/20). 

Unfortunately,  general  real-time  high-resolution  compensation  could  not 
be  demonstrated  with  this  prototype  MSLM,  because  the  parallel -crystal  Fabry- 
Perot  reflections  produce  undesirable11  intensity  modulation  in  an  external 
interferometer.  Now  that  high  quality  prototype  MSLMs  with  oblique  cut  LiNbO 
crystals  are  available,  truly  high  resolution  all-optical  phase  compensation 
should  be  achievable  in  the  near  future. 
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PHASE  COMPENSATED  RESULTS  WITH  ALL-OPTICAL  IPL/MSLM  SYSTEM 


UNCOMPENSATED  COMPENSATED  COMPENSATED 


(A)  (B)  (C) 


exposure:  1/8  sec  1/8  sec  1.5  sec 

(A)  OPEN  LOOP  I NTERFEROGRAM  OF  PHASE  DISTORTION  *0(r)  0F 
ELECTRO-OPTIC  CRYSTAL  IN  AN  EARLY  MSLM. 

(B)  CLOSED  LOOP  (IPL)  INTERFEROGRAM  ILLUSTRATING  PHASE  COMPENSATION 
OF  THE  DISTORTION  SHOWN  IN  (A). 

(C)  LONG-TERM  EXPOSURE  OF  (B)  TO  SHOW  RESIDUAL  ERROR. 

NOTE!  FRINGES  IN  LOWER  LEFT  OF  (A)  AND  (C)  DUE  TO  INTERFERENCE 
BETWEEN  MCP  AND  CRYSTAL  SURFACES  AND  ARE  NOT  WITHIN  THE 
OPTICAL  FEEDBACK  PATH. 

(NO  DIELECTRIC  MIRROR  USED  IN  THIS  DEVICE.) 

Fig  11.  Phase  compensation  results  with  the  all  optical  IPL/MSLM  system. 
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APPENDIX  I 

Optical  information  processing  characteristics  of  the 
microchannel  spatial  light  modulator 

C.  Warde,  A.  M.  Weiss,  A.  D.  Fisher,  and  J,  I.  Thackara 


The  microchannel  spatial  light  modulator  (MSLM)  is  a  versatile,  highly  sensitive,  and  optically  addressed 
modulator  that  is  well  suited  for  low-light  level  real-time  optical  information  processing.  The  image  pro¬ 
cessing  operations  that  can  be  achieved  with  the  MSLM  include  contrast  reversal,  contrast  enhancement, 
edge  enhancement,  image  addition  and  subtraction,  analog  and  digital  intensity  level  thresholding,  and  bi¬ 
nary  level  logic  operations  such  as  AND,  OR,  EXCLUSIVE  OR,  and  NOR.  Several  of  these  operations  are  dem¬ 
onstrated  herein.  Recent  prototype  MSLMs  have  exhibited  a  halfwave  exposure  of  2.2  nj/cm2,  an  optical 
information  storage  time  of  more  than  two  months,  and  a  framing  rate  of  40  Hz  with  full  modulation  depth 
(200  Hz  with  20%  modulation  depth).  The  role  of  secondary  electron  emission  in  the  operation  of  the 
MSLM  is  discussed,  and  design  modifications  that  would  yield  a  spatial  resolution  of  ~10  cycles/mm  at  the 
50%  point  on  an  MTF  curve  are  proposed. 


I.  Introduction 

The  full  exploitation  of  the  inherent  speed  and  par¬ 
allel  processing  capability  of  optical  systems  has  been 
impeded  by  the  lack  of  satisfactory  real-time  reusable 
devices  that  can  be  used  as  substitutes  for  photographic 
film.  Such  devices  are  called  spatial  light  modula¬ 
tors. 

The  microchannel  spatial  light  modulator  (MSLM) 
is  a  relatively  new,  versatile,  highly  sensitive,  and  op¬ 
tically  addressed  device  that  is  well  suited  for  low-light 
level  real-time  optical  information  processing.  In  ad¬ 
dition  to  real-time  operation,  it  can  be  designed  to  ex¬ 
hibit  long  term  (months)  optical  information  storage 
simultaneously.  A  variety  of  sophisticated  image 
processing  operations  can  be  achieved  with  this  device. 
These  include  contrast  reversal,  contrast  enhancement, 
edge  enhancement,  image  addition  and  subtraction, 
analog  and  digital  intensity  level  thresholding,  as  well 
as  binary-level  logic  operations  such  as  and,  OR,  EX¬ 
CLUSIVE  OR,  and  NOR. 

The  MSLM  belongs  to  a  class  of  spatial  light  modu¬ 
lators  that  converts  an  optical  image  into  a  charge  dis- 
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tribution  at  a  crystal-insulating  layer  interface.  The 
charge  distribution  generates  a  voltage  drop  across  the 
crystal  at  the  expense  of  the  voltage  dropped  across  the 
insulating  layer.  The  resulting  electric  field  in  the 
crystal  modulates  the  refractive  index  of  the  crystal. 
The  MSLM  is  unique  in  that  it  employs  a  microchannel 
plate  to  amplify  the  charge  distribution  before  depos¬ 
iting  it  on  the  crystal;  therefore,  the  insulating  layer  in 
the  MSLM  is  constrained  to  be  a  vacuum  gap. 

It  is  the  intent  of  this  paper  to  describe  the  role  of 
secondary  electron  emission  in  the  operation  of  the 
MSLM,  to  give  an  update  of  its  characteristics,  and  to 
describe  and  illustrate  some  of  the  important  image 
processing  operations  that  can  be  achieved  with  this 
device.  Earlier  stages  of  development  of  the  MSLM 
have  been  reported.1'3  Other  real-time  optically  ad¬ 
dressed  spatial  light  modulators  are  reviewed  else¬ 
where.4-6 

II.  Device  Description 

The  MSLM  illustrated  in  Fig.  1  consists  of  a  photo¬ 
cathode  and  a  microchannel  array  plate  (MCP)  in 
proximity  focus  with  a  planar  acceleration  grid  and  an 
electrooptic  plate.  The  electrooptic  plate  generally 
carries  a  high-resistivity  dielectric  mirror  on  one  side 
and  a  transparent  conducting  electrode  on  the  other. 
The  MCP  is  an  array  of  small  (~10-mui)  semiconduct¬ 
ing-glass  pores,  each  functioning  as  a  continuous  dynode 
electron  multiplier.7,8  Spatially  uniform  gains  (±5% 
variation)  of  104  with  a  single  MCP  or  107  with  two 
cascaded  MCPs  (chevron  configuration)  are  typical.9 

In  operation,  the  write  light  (coherent  or  incoherent) 
incident  on  the  photocathode  creates  an  electron  image 
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WRITE  LIGHT 


Fig.  1.  MicroChannel  spatial  light  modulator. 


which  is  amplified  by  the  MCP  and  proximity  focused 
onto  the  surface  of  the  electrooptic  plate.  The  resulting 
surface  charge  distribution  creates  a  spatially  varying 
electric  field  that  modulates  the  refractive  index  of  the 
electrooptic  plate  (crystal).  Thus  the  readout  light, 
which  makes  a  double  pass  through  the  crystal,  is  spa¬ 
tially  phase  modulated. 

The  surface  charge  distribution  may  be  either  nega¬ 
tive  or  positive  depending  on  the  operating  conditions. 
Thus,  two  write  modes  are  possible.  In  the  electron- 
deposition  write  mode,  the  image  is  written  by  adding 
the  appropriate  electron  distribution  to  the  surface  of 
the  crystal;  the  image  may  be  erased  by  secondary 
electron  emission.  In  the  inverted  or  secondary  emis¬ 
sion  write  mode,  the  image  is  written  by  employing 
secondary  electron  emission  to  remove  the  appropriate 
electron  distribution  from  a  uniformly  charged  crystal 
In  this  case,  erasure  can  be  accomplished  by  simply 
adding  electrons  to  the  crystal. 

Optical  information  that  has  been  written  into  the 
MSLM  may  be  read  out  either  interferometrically  with 
schlieren  optics  or  with  crossed  polarizers.  The  possible 
readout  modes  depend  on  the  symmetry  and  cut  of  the 
electrooptic  crystaL  The  readout  mode  that  is  used  in 
the  prototype  devices  employs  the  interference  between 
the  front  and  back  surface  reflections  from  a  plane- 
parallel  electrooptic  crystal.  If  dielectric  mirrors  are 
deposited  on  both  sides  of  such  a  plane-parallel  crystal 
to  form  a  Fabry-Perot  interferometer  of  finesse  F,  a 
complete  off-to-on  transition  can  be  achieved  with  ap¬ 
proximately  1  IF  times  the  halfwave  surface  charge 
density  of  the  crystal. 

The  MSLM  can  be  either  framed  through  a  discrete 
sequence  of  write,  read,  and  erase  periods  or  operated 
in  a  continuous  mode.  In  the  framed  mode,  the  optical 
write  signal  is  integrated  and  stored  in  the  form  of  a 
/  charge  distribution's  removed.  In  the  continuous 
mode  of  operation,  which  can  be  achieved  with  low  re¬ 
sistivity  crystals,  the  instantaneous  modulation  is 
proportional  to  the  write  image  intensity  for  all  tem¬ 
poral  variations  within  the  bandwidth  of  the  system. 

It  is  the  wide  choice  of  photocathodes  and  elec¬ 
trooptic  materials  that  gives  rise  to  the  spectral  versa¬ 
tility  of  the  MSLM  and  permits  its  characteristics  to  be 
tailored  for  specific  applications.  For  example,  the 
wavelength  of  the  write  signal  is  limited  only  by  the 


Fig.  2.  Electrical-circuit  model  of  an  MSLM. 


spectral  response  of  existing  photocathodes.  No  pho¬ 
tocathode  is  needed  in  the  extreme  ultraviolet  and  x-ray 
regions  since  the  MCP  responds  directly  to  this  radia¬ 
tion.1011  The  readout  light  is  limited  by  the  spectral 
transmittance  of  electrooptic  materials  (the  range  be¬ 
tween  the  middle  ultraviolet  and  the  middle  in¬ 
frared). 

Because  of  its  large  electron  gain  and  its  ability  to 
integrate  charge,  the  MSLM  is  capable  of  achieving  the 
high  photon  sensitivities  of  image  intensifier  tubes. 
Devices  employing  one  MCP  are  expected  to  be  able  to 
record  visible  images  with  intensities  as  low  as  10" 14 
W/cm2,  the  primary  limitation  being  MCP  dark  cur¬ 
rent. 

III.  Theory  of  Operation 

An  electrical-circuit  model  of  an  MSLM  is  shown  in 
Fig.  2.  The  electrooptic  crystals  (LiTaOa  and  LiNbOs) 
that  have  been  employed  in  the  prototype  devices  have 
such  high  resistivity  that  Rx  may  be  omitted.  In  this 
case  all  changes  in  Vb  are  divided  capacitively  between 
the  crystal  and  the  gap.  That  is, 


Vb  -  V,  +  V,  +  V., 

(1) 

±V,--^—AVb. 

(2a) 

c,  +  cx 

AVx«— ^-AVV 

(2b) 

C,  +  Cx 

In  general,  the  capacitance  of  the  crystal  is  much  higher 
than  that  of  the  gap,  so  a  greater  fraction  of  the  change 
in  Vb  is  coupled  across  the  gap  than  across  the  crystal. 
The  dependent  current  source  of  density  Jg  in  Fig.  2  is 
the  net  current  density  in  the  gap.  It  is  the  sum  of  the 
primary  current  density  Jp  and  the  secondary  emission 
current  density «/,.  The  write  and  erase  characteristics 
of  the  MSLM  are  determined,  therefore,  by  the  trans¬ 
port,  electrooptical.  and  secondary  emission  properties 
of  the  crystal. 

A.  Monoenergetic  System 

The  classical  case  of  secondary  emission  from  a 
floating  target,  at  potential  Vg,  that  is  addressed  with 
a  monoenergetic  electron  beam  from  a  cathode  at  po¬ 
tential  -V*  (V*  >  0)  is  discussed  at  length  by  Kazan 
and  Knoll.12  The  energy  of  the  primary  electrons  ar¬ 
riving  at  the  surface  of  the  target  of  such  a  system  [see 
Fig.  3(a)]  is  Ep  =  e{Vg  +  Vg ),  where  e  is  the  magnitude 
of  the  electronic  charge. 


After  the  aodulatlon  is  read  out  the  charge  distribution 
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Fig.  3(a).  Idealized,  monoenergetic,  electron-beam  addressed  in¬ 
sulating  target  system  in  which  the  secondary  electrons  are  collected 
by  a  grid  at  ground  potential. 


Fig.  4.  State-space  plot  illustrating  the  behavior  of  a  monoenergetic, 
electron-beam  addressed,  insulating  target  system. 


3(b).  Sketch  of  the  ratio  of  collected  secondary  to  primary  electrons, 
h'{  V*),  for  the  floating  surface  of  an  insulating  target  system.  Positive 
values  of  Vg  are  represented  by  the  solid  lines  and  negative  values  by 
the  broken  lines. 


A  sketch  of  the  ratio  of  collected  secondary  to  primary 
electrons,  b'(  V*),  with  Vg  as  a  parameter,  is  shown  in 
Fig.  3(b).  The  curve  corresponding  to  Vg  =  0  is  the 
conventional  secondary  electron  emission  coefficient 
5(  V*),  and  its  first  and  second  crossover  voltages  are 
labelled  Vt  and  V2.  The  corresponding  primary 
crossover  energies  are  E\  =  eVx  and  £2  =  eV2.  For 
insulators,  the  first  crossover  energy  E\  typically  falls 
in  the  10-100-eo  range  and  the  peak  value  of  5(V*) 
generally  varies  between  1.5  and  25. 12 

The  monoenergetic  state-space  plot  shown  in  Fig.  4 
has  been  adapted  from  Ref.  12.  It  shows  how  the  initial 
kinetic  energy  eV*  of  the  primary  electrons  and  the 
potential  Vg  of  the  floating  surface  must  be  adjusted  to 
add  or  remove  electrons  from  the  target.  External 
control  of  Vg  can  be  achieved  through  V&  in  accordance 
with  Eq.  (2a).  In  Fig.  4,  E\  and  E2  are  the  families  of 
first  and  second  crossover  energies  respectively  of  5'(  V*) 
for  Vg  <  0.  Generally,  E\  is  relatively  independent  of 
Vg  (i.e.,  E\  =  Ei),  but  E\  varies  with  Vg. 12 

When  the  monoenergetic  system  is  operated  in  the 
regions  where  h'  >  1,  the  target  loses  electrons  and  Vg 
decreases,  whereas  when  b'  <  1  the  target  gains  electrons 
and  Vg  increases.  When  5'  =  1,  the  surface  of  the  target 
stabilizes  at  some  equilibrium  voltage  V,,q.  The  plot 
of  Fig.  4  shows  that  there  is  one  line  of  stable  equilibria 
and  one  line  of  unstable  equilibria  that  separate  the 
regions  5'  <  1  and  b'  >  1. 

There  are  three  distinct  operating  regions  for  the 


monoenergetic  system.  These  are  separated  by  vertical 
boundaries  (not  shown)  passing  through  the  conven¬ 
tional  crossover  points  A  and  B  in  Fig.  4: 

(1)  If  e  V*  <  £1,  there  is  no  adjustment  of  Vg  that ' 
will  result  in  the  net  removal  of  electrons  from  the  tar¬ 
get.  This  is  because  at  high  Vg  most  of  the  secondary 
electrons  are  attracted  back  to  the  target  (that  is,  o'  < 
1),  and  for  low  values  of  Vg,  the  low  yield  of  secondary 
electrons  directly  results  in  b'  <  L  Thus,  if  initially  Vg 
>  —V*,  electrons  will  accumulate  on  the  target,  and  Vg 
will  fall  until  the  primary  energy  £p  becomes  zero  (that 
is,  Vg  -*  -  Vk  or  the  potential  drop  across  the  target,  Vx 
—  Vi,  +  V*).  This  is  illustrated  by  the  trajectory  T 1  in 
Fig.  4.  In  the  case  where  the  system  is  operated  with 
Vg  <  -  V*,  no  primary  electrons  can  reach  the  target. 
This  condition  is  called  lockout. 

(2)  If  £1  <  eV*  <  £2  and  the  system  is  operated 
with  large  initial  values  of  Vg,  most  of  the  secondary- 
electrons  will  be  attracted  back  to  the  target.  Vg  will 
fall,  and  b'  will  change  until  it  eventually  becomes  unity 
at  which  point  Vg  stabilizes  at  its  equilibrium  value  Veq. 
Generally  Veq  is  only  a  few  volts.  If,  on  the  other  hand, 
the  system  is  operated  such  that  initially  (Ex  -  e  V\ )  < 
eVg  <  eVeq,  then  since  b'  >  1  in  this  region,  electrons 
will  be  removed  from  the  target,  and  V,  will  rise  toward 
V.  Finally,  if  the  system  is  operated  such  that  —eVk 
<  e  Vg  <  (Ex  —  eVk),  then  since  b'  <  1  in  this  region,  the 
target  will  gain  electrons,  and  Vg  will  fall  until  the  pri¬ 
mary  electron  energy  becomes  zero.  The  various  re¬ 
gimes  of  Vg  for  this  operating  region  are  illustrated  by 
the  trajectory  T2  in  Fig.  4. 

(3)  In  the  case  where  eV*  >  £2,  the  behavior  is 
similar  to  case  2  except  that  the  equilibrium  gap  voltage 
is  negative,  and  its  magnitude  can  be  very'  large,  de¬ 
pending  on  the  value  VV 

B.  MSLM  Multienergetic  System 

The  MSLM  differs  from  the  classical  case  in  that  the 
MCP  is  not  a  monoenergetic  source  of  electrons.  The 
electron  current  from  a  typical  MCP  has  a  continuous 
energy  distribution  n0  ( V* )  ranging  between  zero  and 
eVm,  where  —  Vm  (Vm  >  0)  is  the  MCP  bias  voltage. 
Here,  n0(Vk)dVk  is  the  number  of  MCP  output  elec¬ 
trons  per  unit  area  per  unit  time  with  energy  between 
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Fig.  5.  Sketch  of  the  gap  current  density  J,  in  an  MSLM  as  a  func¬ 
tion  of  gap  voltage  Vg  for  constant  write  light  intensity  Iw  and  MCP 
bias  voltage  -  Vn. 


eVk  and  e<  Vk  +  dVk).  The  MCP  output  current  J0  can 
therefore  be  written  as 

Jo  -  ena(Vk)dVk.  (3) 

In  the  linear  operating  regions  of  the  MCP  and  the 
photocathode,  na(  V*)  scales  linearly  with  the  write  light 
intensicy  Iw.  In  terms  of  the  actual  device  parameters, 
J0  in  this  case  is  given  by 

Jo  *  lwi\eGlhv,  (4) 

where  t)  is  the  quantum  efficiency  of  the  photocathode, 
G  is  the  gain  of  the  MCP,  h  is  Planck’s  constant,  and  u 
is  the  frequency  of  the  write  light.  The  primary  current 
Jp  is  defined  at  the  surface  of  the  crystal,  and  in  a 
gridless  device  Jp  =  J„  when  none  of  the  MCP  electrons 
is  locked  out  (repelled  from  the  crystal). 

It  has  been  shown1314  that  n0(Vk)  is  weighted  heavily 
in  favor  of  the  low  energy  electrons.  Typically,  n0  ( Vk ) 
spans  region  I  and  part  of  region  II  in  Fig.  4.  Thus,  in 
the  gridless  device,  the  equilibrium  voltage  is  often 
negative. 

In  devices  employing  an  acceleration  grid  at  a  po¬ 
tential  Va ,  the  primary  current  Jp  is  always  less  than 
the  MCP  output  current  J0.  If  n(Iu.,Vk,Va,Vg )  is  the 
primary  electron  energy  distribution  at  the  surface  of 
the  crystal,  then 

JP  -  j‘oVmen(Iw,Vk,V„,Vg)dVk  (5) 

Additionally,  an  effective  collected  secondary  electron 
emission  ratio  8'  may  also  be  defined  by 

V(Vm,Vo,V,)  •  CVme6'(Vl„Vo,Vg)r,aw,Vk.V<„V,)dVk.  (6) 
Jp  *Jo 

The  gap  current  density  Jg  determines  the  rate  at 
which  electrons  are  added  to  the  crystal.  It  is  given 
by 

Jg  **  Jp{  1  —  6’).  (7) 

Figure  5  is  a  sketch  of  the  gap  current  density  Jg  as  a 
function  of  Vg,  with  Iw,  Va  and  Vm  constant,  for  a  de¬ 
vice  employing  a  grid  and  operated  such  that  Veq  is 
positive. 

The  MSLM  is  operated  so  that  it  is  either  at  the 
equilibrium  point  E  in  Fig.  5  or  it  is  moving  along  one 
of  the  trajectories  of  Fig.  5  toward  equilibrium.  For 
materials  exhibiting  the  linear  electrooptic  effect,  two 
write  modes  are  therefore  possible.  In  the  electron 


deposition  write  mode,  V*  initially  divides  capacitively 
between  the  crystal  and  the  gap,  thereby  raising  Vg  into 
the  region  for  which  8’  <  1 .  Then  when  the  write  light 
is  turned  on,  electrons  are  deposited  on  the  crystal,  Vg 
begins  to  fall  toward  its  equilibrium  value  V’?(J ,  and  the 
crystal  voltage  Vx  rises.  Electron  accumulation  may 
be  terminated  when  the  desired  Vx  is  reached  by  turn¬ 
ing  off  the  write  light.  The  crystal  may  be  erased  via 
secondary  electron  emission  by  flooding  the  photo¬ 
cathode  with  light  an_d  lowering  Vg  (by  lowering  Vb) 
into  the  region  where  8'  >  1.  In  practice,  Vb  is  ramped 
downward  no  faster  than  Vg  can  be  restored  to  Veq  by 
electron  removal,  so  as  to  avoid  lockout. 

In  the  secondary  electron  emission  write  mode  with 
the  write  light  on,  Vb  is  ramped  downward  at  a  rate  Vb 
chosen  so  that  lockout  is  avoided  and  the  appropriate 
electron  distribution  is  removed  from  the  crystal.  The 
secondary  emission  write  process  terminates  when  ei¬ 
ther  the  write  light  is  turned  off,  or  Vb  reaches  its  final 
value  and  Vg  returns  to  its  equilibrium  value.  To  erase 
the  information,  Vb  can  be  increased  abruptly  and  the 
photocathode  flooded  with  light  so  that  the  crystal  gains 
electrons. 

During  the  removal  of  electrons  by  secondary  emis¬ 
sion  from  the  crystal  in  a  gridless  device,  the  low  energy 
MCP  output  electrons  for  which  5'  <  1  tend  to  increase 
rather  than  decrease  the  electron  density  at  the  crystal 
surface  thereby  resulting  in  inefficient  erasure.  How¬ 
ever,  this  problem  is  overcome  in  devices  employing  an 
acceleration  grid.  The  grid  potential  Va  is  adjusted  to 
shift  the  center  of  gravity  of  n(Iw,Vk,Va,Vg)  so  that  8’ 
is  maximized.  A  measurement  of  the  secondary  elec¬ 
tron  coefficient  for  LiNbC>3  is  reported  in  Ref.  15. 


IV.  Write  and  Erase  Characteristics 


From  the  foregoing  analysis,  it  is  clear  that  for  both 
write  modes  the  fastest  framing  rate  is  achieved  by  op¬ 
erating  the  device  during  the  electron  deposition  period 
with  Vg  large  enough  to  have  Jg  cs  J0  (see  Fig.  5)  and 
by  removing  electrons  with  the  device  operated  at  the 
point  M  where  the  crystal  loses  electrons  at  the  fastest 
rate.  Thus,  the  minimum  halfwave  electron  deposition 
time  td*  for  a  high-resistivity  crystal  is  given  by 


tdw 


_  VrU(o/lx 
Jo  Jo 


*81 


and  the  minimum  halfwave  electron  removal  time  trir 
is  given  by 


tn  -  y&m  -  l)-‘.  (9) 

where  8m  is  the  maximum  value  of  8\  f0  is  the  permit¬ 
tivity  of  free  space.  and  lx  are  the  dielectric  constant 
and  the  thickness  of  the  crystal  parallel  to  the  charge- 
induced  electric  field,  respectively,  and  <xT  and  VT  are 
the  reflection  halfwave  surface  charge  density  and 
halfwave  voltage,  respectively,  of  the  crystal. 

The  reflection  halfwave  surface  charge  density  <rT  is 
the  charge  density  at  the  surface  of  the  electrooptic 
crystal  that  produces  a  phase  shift  of  7r  rad  in  the 
readout  beam  for  a  double  pass  through  the  crystal. 
For  r-cut  LiNbC>3  and  LiTaC>3  operated  in  the  longi¬ 
tudinal  mode,  the  electrically  induced  phase  change  is 
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independent  of  the  polarization  of  the  light  (i.e.,  A‘l>  = 
A#,  =  A4»y),  and 

a,  “  «  i<oX/2no''i3/, ,  (10) 

where  X  is  the  wavelength  of  the  readout  light,  n0  is  the 
ordinary  refractive  index  of  the  crystal,  and  ri3  the 
electrooptic  coefficient  of  interest.  For  a  1-mm  thick 
LiNb03  crystal,  <r,  is  approximately  0.8  X  10-7  C/cm2 
at  6328  A. 

The  write  and  erase  times  are  generally  limited  by  the 
maximum  MCP  current  available. 


V.  Exposure  Sensitivity 

The  halfwave  exposure  S„  of  the  device  is  the  write 
light  energy  per  unit  area  that  is  required  to  change  the 
phase  of  the  readout  light  by  ir  rad  and  is  given,  for  both 
write  modes,  by 

S„  -  £’lu(t)dt.  (11) 

where  tr  is  given  by  either  Eq.  (8)  or  (9). 


VI.  Spatial  Resolution 

The  spatial  resolution  of  the  MSLM  is  determined 
by  the  fidelity  of  the  proximity  focusing16  and  by  the 
fringing  of  the  electric  fields  in  the  electrooptic  crystal. 
To  maintain  good  proximity  focusing,  the  longitudinal 
velocity  component  of  the  primary  electrons  must  be 
much  greater  than  the  transverse  component.  This 
condition  is  achieved  in  the  gridless  device  by  setting 
Vb  a  few  hundred  volts  above  the  desired  voltage  swing 
of  the  crystal.  In  devices  employing  an  acceleration 
grid,  the  grid  voltage  V„ ,  which  optimizes  electron  re¬ 
moval  from  the  crystal,  is  generally  enough  to  impart 
sufficient  longitudinal  kinetic  energy  to  the  primary 
electrons  to  maintain  good  proximity  focusing. 

The  case  of  crystal-limited  resolution,  which  applies 
to  all  MSLMs  built  so  far,  has  been  discussed  at  length 
by  Roach17  and  generalized  by  Owechko  and  Tanguay.18 
For  structures  similar  to  the  MSLM,  Roach  has  shown 
that  the  crystal-limited  resolution,  as  measured  by  the 
spatial  bandwidth  w i  (MTF  roll-over  frequency)  is 
given  by 


.  tiH* +  u!ii  , 
»IUj«x)‘/2  +  u\ 


(12) 


where  lg  is  the  gap  width,  f  ±  is  the  dielectric  constant 
perpendicular  to  the  direction  of  the  charge-induced 
electric  field,  tg  is  the  dielectric  constant  of  the  vacuum 
gap,  and  eg  ana  lx  are  as  defined  in  Eq.  (8). 

From  Eq.  (12)  it  is  clear  that,  in  principle,  the  spatial 
bandwidth  can  be  increased  by  reducing  either  the  gap 
width  lg  or  the  crystal  thickness  lx.  However,  the 
product  of  the  crystal-limited  spatial  bandwidth  w\  and 
the  low-frequency  halfwave  charge  sensitivity  S„(0)  of 
such  devices  is  a  constant  which  does  not  depend  on 
device  dimensions. 17  Thus,  any  improvement  in  spatial 
bandwidth  will  be  at  the  expense  of  charge  sensitivity. 
For  the  MSLM,  SJO)  is  given  by 


S,(0) 


(13) 


For  convenience  all  MSLMs  built  so  far  have  em¬ 
ployed  thick  (0.5-mm)  crystals  and  wide  (~0.5-mm) 


0.1  1.0  10  100 
w(cycles/mm) 

Fig.  6.  Calculated  normalized  charge  sensitivity  S„(ivl/S„(0)  vs 
spatial  frequency  w  for  an  oblique-cut  LiNbO.-j  MSLM  with  crystal 
thickness  and  gap  width  as  parameters.  The  crystal  cut  is  such  that 
the  normal  to  the  surface  makes  an  angle  of  55°  with  the  z  axis  and 
lies  in  the  plane  containing  the  z  and  —y  axes. 


gaps.  Future  LiNb03  and  LiTa03  devices  will  employ 
crystals  75-250  gm  thick  and  gaps  of  ~100-500  jam.  A 
z  -cut  LiTa03  device  with  a  75-jam  thick  crystal  and  a 
100-Mm  gap  will  have  an  MTF  rollover  of  about  4.7  cy¬ 
cles/mm  which  corresponds  to  a  spatial  resolution  of 
better  than  10  cycles/mm  at  the  50%  point  on  an  MTF 
curve.  The  minimum  crystal  thickness  is  limited  by 
dielectric  breakdown  and  by  the  difficulties  associated 
with  polishing  thin  crystals;  the  minimum  gap  width 
must  be  large  enough  to  prevent  arcing  during  opera¬ 
tion. 

The  MSLM  can  be  further  optimized  for  both  in¬ 
tensity  and  phase  modulation  by  employing  a  crystal 
cut  that  minimizes  the  halfwave  surface  charge  density. 
For  LiNb03,  this  cut  is  such  that  the  normal  to  the 
crystal  surface  lies  in  the  plane  containing  the  z  and  —y 
axes  and  makes  an  angle  of  55°  with  the  z  axis.19  With 
this  cut,  the  electrically  induced  birefringence  alters  the 
polarization  of  the  readout  light,  and  the  reflection 
halfwave  voltage  falls  to  a  minimum  value  of  ~1250  V. 
Figure  6  shows  how  the  calculated  spatial  resolution  of 
such  a  device  depends  on  crystal  thickness  and  gap 
width. 

VII.  Storage 

Stored  charge  images  will  diffuse  with  a  characteristic 
relaxation  time  that  depends  on  the  finite  surface  and 
bulk  conductivities  of  the  crystal  or  the  dielectric  mir¬ 
ror.  The  device,  therefore,  must  be  framed  through  the 
write,  read,  and  erase  periods  in  a  time  that  is  much 
shorter  than  the  charge  relaxation  time. 

VIII.  Prototype  Devices 

Two  simple  prototype  devices  have  been  built.  The 
first  is  a  vacuum  demountable  device  employing  a 
25-mm  diam  MCP  in  proximity  focus  with  a  nominally 
plane-parallel  0.5-mm  thick  crystal  that  is  coated  on 
one  side  with  a  transparent  electrode  of  indium  tin 
oxide.  Three  z  -cut  crystals  have  been  evaluated  in  this 


2070  APPLIED  OPTICS  /  Vol.  20.  No.  12  /  15  June  1981 


WRITE 

LIGHT 


—  Cryttol  L'tcroat 
MCP  out 


Fig.  7.  Vacuum  sealed,  visible  photocathode,  prototype  MSLM 
employing  electron  optics  between  the  photocathode  and  the  MCP 
and  proximity  focusing  between  the  MCP  and  the  crystal. 
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Fig.  9.  Interferometric  readout  intensity  characteristic  for  an  un¬ 
coated  LiTa03  or  LiNbCb  crystal  (ft  *  0.15)  and  for  a  Fabry-Perot 
(R  =  0.80)  crystal. 


Fig.  8.  Setup,  with  interferometric  readout,  used  to  measure  the 
write  and  erase  times,  the  spatial  resolution,  the  image  storage  char¬ 
acteristics.  and  the  operations  achievable  with  the  MSLM. 


device:  One  LiTa03  and  one  LiNb03  crystal  without 
dielectric  mirrors,  and  one  LiNbOg  crystal  with  di¬ 
electric  mirrors  on  both  sides  (Fabry  Perot  configura¬ 
tion).  In  each  case,  the  jap  between  the  MCP  and  the 
crystal  was  0.5  mm  wide.  A  planar  acceleration  grid  was 
sometimes  inserted  into  the  gap.  The  Inconel  electrode 
on  the  input  face  of  the  MCP  served  as  the  photocath¬ 
ode.  This  device  was  activated  with  ultraviolet  write 
light  (254  nm)  from  a  mercury  lamp. 

The  second  device  is  vacuum  sealed  and  employs  a 
visible  bialkali  photocathode  with  6%  quantum  effi¬ 
ciency  at  655-nm  wavelength.  The  gap  in  this  device 
is  0.5  mm  wide,  and  the  electrooptic  plate  is  a  nominally 
plane-parallel  0.5-mm  thick  2 -cut  LiTa03  crystal 
without  a  dielectic  mirror.  There  is  no  grid  in  this  de¬ 
vice.  To  simplify  fabrication,  this  device  was  built  by 
modifying  a  Hamamatsu  Corp.  image  intensifier  tube. 
In  this  design,  electron  optics  are  used  to  image  elec¬ 
trons  from  the  photocathode  onto  the  MCP,  and  prox¬ 
imity  focusing  is  employed  between  the  MCP  and  the 
crystal.  The  vacuum  sealed,  visible  photocathode  de¬ 
vice  is  shown  schematically  in  Fig.  7.  An  LED  was  used 
as  the  source  of  write  light  for  this  device. 

The  arrangement  shown  in  Fig.  8  was  used  to  measure 
the  write  and  erase  times,  the  spatial  resolution,  the 
image  storage  characteristics,  and  the  image  processing 
operations  achievable  with  these  devices.  An  inter¬ 
ferometric  readout  scheme  employing  the  Haidinger 


interference  patterns  produced  by  multiple  reflections 
between  the  front  and  back  surfaces  of  the  plane-par¬ 
allel  crystals  was  used.  The  surface  reflectivity  of  un¬ 
coated  LiTa03  and  LiNb03  (0.15  at  normal  incidence) 
provided  sufficient  readout  image  contrast  for  experi¬ 
mental  purposes  (see  Fig.  9). 

The  demountable  LiTa03  device  without  a  dielectric 
mirror  and  with  a  grid  in  the  gap  was  cycled  at  20 
frames/sec  with  full  modulation  depth  and  at  100 
frames/sec  with  20%  modulation  depth.  The  de¬ 
mountable  Fabry  Perot  LiNb03  device,  employing  an 
acceleration  grid,  was  cycled  at  40  frames/sec  with  full 
modulation  depth,  and  up  to  200  frames/sec  with  20% 
modulation  depth.  The  MCP  in  this  device  had  a  strip 
current  of  20  iiA  (at  Vm  =  1  kV).  MCP  manufacturers 
are  now  developing  MCPs  with  strip  currents  in  excess 
of  100  /uA  (at  Vm  -  1  kV),  so  it  is  expected  that  the 
speed  of  the  device  can  be  increased  by  at  least  a  factor 
of  5  in  the  near  future. 

The  measured  halfwave  exposure  of  the  visible  pho¬ 
tocathode  device  for  the  electron  deposition  write  mode 
was  2.2  nj/cm2  at  655-nm  wavelength,  making  it  ~3 
orders  of  magnitude  more  sensitive  than  the  other 
real-time,  optically  addressed  spatial  light  modula¬ 
tors.4-6  The  gain  of  the  MCP  in  this  device  was  less 
than  1000.  Therefore,  with  the  standard  MCP  gain  of 
104,  the  halfwave  exposure  of  this  device  should  be 
comparable  to  the  calculated  value  of  0.3  nJ/cm2.1 

With  the  visible  photocathode  LiTa03  device,  optical 
information  storage  times  of  up  to  two  months  with 
negligible  charge  diffusion  have  been  observed.  Be¬ 
cause  of  the  thick  (0.5-mm)  crystals  and  the  wide 
(~0.5-mm)  gaps  used  in  the  prototype  devices,  their 
spatial  resolution  was  only  ~2  cycles/mm  at  the  50% 
point  on  an  MTF  curve.  However,  as  shown  in  Sec.  VI, 
the  resolution  can  be  improved  to  ~10  cycles/mm  at  the 
50%  point  by  employing  75  to  100  nm  thick  crystals. 

IX.  Achieveable  Image  Processing  Operations 

By  exploiting  the  secondary  electron  emission  char¬ 
acteristics  of  the  surface  that  stores  the  charge  distri¬ 
bution,  a  variety  of  sophisticated  image  processing  op¬ 
erations  are  possible.  These  include  analog  and  digital 
intensity  thresholding  on  stored  images,  edge  en¬ 
hancement,  contrast  enhancement,  contrast  reversal. 
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10 (A), (B)  CONTRAST  REVERSAL 


11 (A), (B)  IMAGE  ADDITION 


12 (A)-(C)  IMAGE  SUBTRACTION 
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13 (AMC)  ANALOG  AND  DIGITAL  THRESHOLDING 
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15 (AMB)  m  OPERATION 

Pigs.  10-15.  Some  of  the  image  processing  operations  achieved  with 
the  visible  photocathode  prototype  MSLM. 

image  addition  and  subtraction,  and  binary-level  logic 
operations  such  as  and,  OR,  exclusive  or.  and  nor. 
The  binary  (two-level  irradiance)  logic  operations  are 
best  performed  in  a  device  employing  a  grid.  The  grid 
accelerates  the  MCP  electrons  and  sets  the  maximum 
charge  density  level.  Thus,  any  point  on  the  crystal  can 
be  written  to  the  full  charge  density  without  a  loss  of 
proximity  focusing. 

The  MSLM  interferometric  readout  intensity  char¬ 
acteristic  Ir  as  a  function  of  crystal  surface  charge 
density  a  is  shown  in  Fig.  9.  The  solid  curve  is  for  an 
uncoated  LiTaCh  or  LiNb(>3  crystal  such  as  that  used 
in  the  visible  photocathode,  prototype  MSLM.  The 


reflectivity  R  of  the  uncoated  surfaces  of  these  crystals 
is  ~0.15.  The  dashed  curve  is  typical  of  a  Fabry- Perot 
crystal  with  equal  surface  reflectivities.  The  point  Po 
is  a  null  reflection  point,  and  <to  is  the  surface  charge 
density  required  to  bias  the  crystal  at  Po-  The  points 
H+  and  H—  are  half-power  points,  and  the  points  Pi_ 
and  Pi+  are  the  90%  of  peak  transmission  points.  The 
point  Pr  is  a  peak  transmission  point  with  a  corre¬ 
sponding  crystal  surface  charge  density  of  <ro  +  ov. 
Similar  points  can  be  identified  on  the  R  =  0.15  char¬ 
acteristic.  Notice  the  improved  charge  sensitivity  for 
the  off-to-on  transition  for  a  Fabry-Perot  device. 

All  the  image  processing  operations  illustrated  here 
were  achieved  with  the  vacuum  sealed,  visible  photo¬ 
cathode,  prototype  MSLM  using  the  arrangement 
shown  in  Fig.  8.  The  fine  fringe  pattern  that  modulates 
all  the  readout  images  in  Figs.  10-15  is  a  result  of  in¬ 
terference  caused  by  reflections  from  the  contiguous 
surfaces  of  the  MCP  and  the  crystal.  This  pattern  will 
disappear  when  crystals  with  dielectric  mirrors  are 
used. 

A.  Contrast  Reversal 

As  was  pointed  out  earlier,  images  can  be  written  into 
the  MSLM  as  either  positive  or  negative  charge  distri¬ 
butions.  And  since  both  LiNbC>3  and  LiTaC>3  exhibit 
the  linear  electrooptic  effect,  contrast  reversed  images 
can  be  written  by  reversing  the  polarity  of  the  surface 
charge  distribution.  This  is  illustrated  in  Figs.  10(A) 
and  (B)  for  a  square  input  optical  image  written  as  a 
negative  and  a  positive  charge  distribution,  respectively, 
in  the  visible  photocathode  MSLM. 

The  operation  of  contrast  reversal  can,  in  principle, 
also  be  performed  on  stored  images.  If,  for  example,  an 
image  is  stored  along  the  path  P0P1+,  contrast  reversal 
can  be  achieved  either  by  adding  a  uniform  electron 
distribution  of  <xH_  —  <jh+  or  by  removing  a  uniform 
electron  density  of  oh+  ~  <?h~  In  practice  it  is  easier 
to  add  than  to  remove  large  quantities  of  charge  uni¬ 
formly. 

B.  Image  Addition  and  Subtraction 

The  addition  of  two  images  is  best  achieved  by  suc¬ 
cessively  recording  both  images  by  the  electron-depo¬ 
sition  mode.  Image  addition  with  the  visible  photo¬ 
cathode  MSLM  is  illustrated  in  Figs.  11(A)  and  (B)  for 
a  set  of  horizontal  line  sources  [image  Ai]  and  a  set  of 
vertical  line  sources  [image  A2]. 

To  implement  image  subtraction,  the  base  line 
readout  intensity  of  the  MSLM  is  first  set  at  a  half¬ 
power  point.  This  offset  is  necessary  to  ensure  that  all 
points  of  the  result  (after  subtraction)  are  positive.  If, 
for  example,  the  base  line  half-power  point  H+  is  used, 
the  first  image  Aj  is  written  by  secondary  electron 
emission,  and  the  second  image  A2  is  written  by  electron 
deposition.  The  exposures  of  Ai  and  A2  must  be  con¬ 
trolled  so  that  all  parts  of  both  images  are  written  within 
the  most  linear  part  of  the  readout  intensity  charac¬ 
teristic.  The  operations  leading  to  A2  -  A(  are  illus¬ 
trated  in  Figs.  12(  A)-(C)  for  a  horizontal  bar  image  At 
which  intersects  a  vertical  bar  image  A2. 
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C.  Analog  Thresholding 

Consider  two  areas  A'  and  A"  within  an  electron 
distribution  on  the  surface  of  the  crystal  with  electron 
densities  o'  and  o",  respectively.  If  a”  is  sufficiently 
greater  than  o ',  then  Vb  can  be  adjusted  so  that,  when 
the  photocathode  is  uniformly  illuminated,  primary 
electrons  from  the  MCP  are  repelled  from  A"  but  not 
from  A'.  Under  these  conditions  A'  and  all  areas  with 
electron  density  less  than  o'  will  be  erased,  while  A"  and 
all  other  areas  for  which  o>  a”  will  be  unaffected.  This 
operation  (analog  level  thresholding)  corresponds  to 
erasing  all  the  parts  of  an  image  for  which  the  intensity 
falls  below  some  threshold  intensity  It.  An  example  of 
analog  thresholding  is  illustrated  in  Figs.  13(A)  and 
(B). 

Here  analog  thresholding  is  used  to  erase  the  lowest 
intensity  level  in  a  three-level  image  leaving  the  higher 
intensity  levels  unaffected. 

Repeated  applications  of  this  operation  (multiple- 
level  thresholding)  can  be  used  to  analyze  the  spatial 
distribution  of  intensities  in  an  image.  Analog 
thresholding  is  useful  in  the  implementation  of  contrast 
enhancement  and  the  logic  operations  AND  and  XOR. 

O.  Stored  Image  Digital  Thresholding  (Hard  Clipping) 

In  binary-level  optical  logic  operations  it  is  often 
desirable  to  digitize  images  by  means  of  an  operation  in 
which  all  intensities  below  some  chosen  threshold  level 
It  are  transformed  to  the  same  lower  level  )e.g.,  S  =  0 
in  Fig.  9).  and  all  those  above  It  are  restored  to  the  same 
upper  level  (e.g..  S  *  1).  This  operation,  also  called 
hard  clipping,  can  be  achieved  in  a  three-step  process 
of  which  the  first  step  is  analog  thresholding.  In  the 
second  step,  the  photocathode  is  flooded  uniformly  with 
light  and  all  the  remaining  nonzero  intensities  are 
transformed  to  the  same  saturation  intensity  level  set 
by  the  grid  voltage.  Then  analog  thresholding  is  re¬ 
peated  to  restore  the  lower  level  intensities  to  the  S  = 
0  state.  Figure  13(C)  illustrates  the  hard  clipping  result 
derived  after  the  analog  thresholding  operation  [Figs. 
13(A)  and  (B)]  on  the  three-level  system. 

E.  Real-Time  Hard  Clipping 

This  operation  is  most  conveniently  performed  in  the 
secondary  emission  write  mode.  The  MSLM  is  initially 
biased  in  the  OFF  state  with  charge  density  oo-  Then 
Vb  is  ramped  downward  (from  higher  to  lower  voltage) 
at  a  rate  Vb  with  the  optical  image  incident  on  the 
photocathode.  All  the  intensities  below  some  threshold 
level  I,  (which  depends  on  Vb)  will  be  barely  recorded 
because,  at  their  locations,  electrons  cannot  be  removed 
fast  enough  from  the  crystal  to  prevent  the  gap  field 
from  eventually  repelling  all  future  primary  electrons. 
Thus,  all  intensities  below  threshold  will  remain  at  a  low 
intensity  level  (S  *  0  state).  Because  the  slope  of  the 
Vx  vs  Iw  characteristic  is  known20  to  be  very  steep  be¬ 
tween  threshold  and  saturation,  virtually  all  write  light 
intensities  above  threshold  are  written  to  the  maximum 
charge  density  which  is  set  to  coincide  with  the  S  *  1 
state. 


This  operation  can  lead  to  real-time  contrast  en¬ 
hancement  if  the  crystal  has  a  low  halfwave  voltage  (e.g., 
DKDP  at  -51°C,  V,  **  200  V)  and  if  Vb  is  adjusted  so 
that  the  spatial  intensity  variations  of  the  image  are 
mapped  onto  the  knee  of  the  Vx  vs  Iw  curve. 

F.  Contrast  Enhancement 

Contrast  enhancement  of  stored  images  is  a  direct 
consequence  of  analog  thresholding  since  in  this  oper¬ 
ation  the  low-level  intensities  in  an  image  can  be  re¬ 
duced  to  an  even  lower  level,  leaving  the  high-level  in¬ 
tensities  unaffected.  An  example  of  contrast  en¬ 
hancement  using  level  thresholding  is  illustrated  in  Figs. 
14( A)— (C)  with  the  visible  photocathode  MSLM.  In 
Fig.  14(B)  the  stored  image  shown  in  Fig.  14(A)  is 
washed  out  by  the  addition  of  a  uniform  charge  density 
to  the  surface  of  the  crystal.  In  Fig.  14(C)  the  contrast 
is  enhanced  by  the  level  thresholding  operation. 

G.  Stored-lmage  Edge  Enhancement 

Suppose  that  a  stored  image  of  relatively  uniform 
intensity  was  written  with  an  average  charge  density 
level  of  o  on  a  uniform  background  of  charge  density  rr0- 
Edge  enhancement  can  be  performed  on  this  image  by 
flooding  the  photocathode  uniformly  with  light  until  the 
average  charge  density  of  the  image  reaches  a  level  of 
Oq  +  o-it,  (OFF  state),  and  then  using  analog  thresholding 
to  restore  the  background  to  the  level  (To  (OFF  state). 
Under  these  conditions,  the  charge  density  will  fall 
abruptly  from  an  average  value  of  o0  +  <r2)r  inside  to  a 
value  of  <r0  outside  the  image.  Thus,  on  the  average  the 
charge  density  at  the  edges  of  the  image  will  be  at  a  level 
of  (To  +  <rr  which  corresponds  to  the  ON  state.  That  is, 
the  edges  of  the  image  will  be  enhanced. 

H.  Real-Time  Edge  Enhancement 

With  the  modulator  biased  in  the  OFF  state  at  P0  with 
a  uniform  charge  density  oq,  real-time  edge  enhance¬ 
ment  of  an  image  of  relatively  uniform  intensity  can  be 
achieved  by  writing  the  image,  in  the  electron  deposition 
mode,  to  an  average  level  of  <r0  +  <r2jr  (OFF  state).  Again, 
the  edges  of  the  image  will  have  an  average  charge 
density  of  oq  +  or  (ON  state)  and,  therefore,  will  be  en¬ 
hanced  relative  to  the  rest  of  the  image.  Notice  that  the 
modulator  must  be  cycled  through  a  full  wave  of  phase 
modulation  in  order  to  implement  this  operation. 
Because  of  possible  damage  to  the  visible-photocathode 
prototype  modulator  we  have  not  yet  attempted  this 
operation. 

I.  or  Operation 

The  implementation  of  this  operation  depends  on  the 
type  of  device  being  used.  Consider  first  a  Fabry-Perot 
device  with  the  high -reflectivity  readout  characteristic 
shown  in  Fig.  9.  The  base  line  charge  density  is  first  set 
at  <t0  which  corresponds  to  the  S  «  0  state.  If  two  bi- 
nary-level  images  Bi  and  B2  are  now  each  written  one 
at  a  time  with  a  charge  density  <ti,  the  charge  density  at 
the  intersection  of  images  Bj  and  B2  will  be  cr0  +  2oi 
(the  S  »  1  state),  whereas  all  the  other  parts  of  these 
images  will  have  a  charge  density  of  o0  +  o i  (also  the  S 
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=  1  state).  The  background  will  remain  in  the  S  -  0 
state  with  a  charge  density  of  <xo.  The  result  will  be  the 
OR  operation. 

If  the  device  exhibits  the  low-reflectivity  readout 
characteristic  shown  in  Fig.  9,  after  images  Bi  and  Bo 
are  each  written  with  a  charge  density  of  <rj 2,  digital 
thresholding  must  be  performed  to  restore  all  parts  of 
the  composite  image  to  a  charge  density  of  <to  +  0* 
which  corresponds  to  the  S  =  1  state  and  all  parts  of  the 
background  to  a  charge  density  of  a0  which  corresponds 
to  the  S  =  0  state. 

J.  and  Operation 

This  operation  is  implemented  in  the  same  way  for 
both  types  of  device.  In  each  case  the  base  line  charge 
density  is  set  at  cro  (S  =  0  state),  and  the  two  binary- 
level  images  Bi  and  B2  are  written  one  at  a  time  with  a 
charge  density  of  ajP.  Analog  thresholding  is  then 
used  to  remove  all  elecirv  n  densities  that  are  less  than 
or  equal  to  <ro  +  uj2 ;  the  result  will  be  the  AND  opera¬ 
tion  of  Bi  and  B2.  The  AND  operation  is  illustrated  in 
Figs.  15(A)  and  (B)  with  the  visible  photocathode 
MSLM  for  a  horizontal  pair  of  line  sources  [image  Bi] 
and  a  vertical  pair  of  line  sources  [image  Bo]. 

K.  exclusive  or  (xor)  Operation 

This  operation  is  implemented  in  the  same  way  for 
both  types  of  device.  Starting  with  the  modulator  at 
Po  in  the  S  =  0  state  with  a  uniform  charge  density  <r0, 
images  Bi  and  Bo  are  written,  one  at  a  time,  with  a 
charge  density  of  aT.  The  intersection  of  the  two  im¬ 
ages  will  then  be  in  the  S  *  0  state  at  the  point  P0  with 
a  charge  density  of  <ro  +  02*.  whereas  all  the  other  parts 
of  images  Bx  and  B2  will  remain  in  the  S  =  1  state  with 
a  charge  density  of  <to  +  0*.  The  result  will  be  the  XOR 
operation. 

L.  nor  Operation 

With  the  grid  voltage  set  at  V0  +  Vj*.  the  steps 
leading  to  the  XOR  operation  are  performed  first.  The 
photocathode  is  then  uniformly  flooded  with  light  until 
all  the  parts  of  the  composite  image  [Bi  and  Bo]  reach 
the  maximum  charge  density  of  <r0  +  cr2»-  which  corre¬ 
sponds  to  the  S  =  0  state.  The  background  will  then 
have  a  charge  density  of  00  +  0*  (5  *  1  state).  Thus  the 
NOR  operation  is  achieved. 

X.  Comments 

The  spectral  versatility,  the  high  exposure  sensitivity, 
the  long  storage  times,  and  the  wide  variety  of  image 
processing  operations  that  can  be  achieved  by  exploiting 
secondary  electron  emission  are  the  most  outstanding 
features  of  the  MSLM.  Since  none  of  the  characteris¬ 
tics  of  the  present  prototype  devices  is  at  its  funda¬ 
mental  limit,  it  should  be  possible  to  improve  all  aspects 
of  the  performance  of  the  MSLM. 

Specifically,  it  should  be  possible  to  realize  an  expo¬ 
sure  sensitivity  of  better  than  1  nJ/cm-  in  either  a 
LiTaO.i  or  a  LiNbOs  device  by  employing  an  MCP  with 
the  standard  gain  of  104.  Also,  the  present  framing  rate 
of  40  Hz  (full  modulation  depth)  can  be  improved  even 
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further  by  (1)  using  a  higher  strip  current  micr-'channel 
plate,  and  (2)  coating  the  surface  of  the  crystal  with  a 
material  (such  as  MgO)  that  has  a  very  high  secondary 
electron  emission  coefficient.  Additionally,  by  reducing 
the  crystal  thickness  to  ~75  urn,  the  spatial  resolution 
at  the  50%  point  on  an  MTF  curve  should  be  ~10  cv- 
cles/mm.  For  such  thin  crystals  it  will  be  necessary  to 
employ  the  optimum  crystal  cut  so  that  the  operating 
voltages  and.  therefore,  the  probability  of  dielectric 
breakdown  in  the  crystal  can  be  significantly  re¬ 
duced. 

Despite  the  low  resolution  of  the  visible  photocathode 
prototype  MSLM,  it  has  been  demonstrated  that  the 
MSLM  is  also  capable  of  performing  a  wide  variety  of 
important  signal  processing  operations  in  real  time  as 
well  as  on  stored  images.  The  accuracy,  dynamic  range, 
and  linearity  of  these  operations  are  currently  being 
investigated.  Alternative  means  of  implementing  logic 
operations  with  the  MSLM  are  described  in  Ref.  20. 
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Significant  improvements  in  the  spatial  resolution,  framing  speed,  optical  quality,  and  space-domain  image-pro- 
cessing  capabilities  of  the  microchannel  spatial  light  modulator  have  been  realized  by  employing  oblique-cut  rather 
than  z-cut  LiNbOj  crystals,  high-strip-current  <  250-pA)  microchannel  plates,  and  an  acceleration  grid  in  the  gap 
of  the  device.  In  particular,  a  prototype  device  employing  a  330-Mm-thick.  optimum-cut  (rotated  35°  from  the  z 
axial  LiNbOi  crystal  exhibited  a  framing  rate  in  excess  of  30  Hz  with  full  modulation  depth  and  a  spatial  resolution 
of  ~1.9  cycles/ mm  at  50%  contrast  i  ~  10  cycles/  mm  at  10%  contrast  1.  Additionally,  four-level  to  two-level  intensity 
image  conversion,  contrast  reversal,  contrast  enhancement,  edge  enhancement,  and  the  binary-level  operations 
AND,  NAN’D,  OR,  NOR.  XOR.  and  NXOR  have  been  demonstrated  by  operating  the  device  in  its  space-domain 
image-processing  mode. 


The  microchannel  spatial  light  modulator  (MSLM)  is 
a  versatile,  retd-time  optical  signal-  and  image-pro¬ 
cessing  device  that  exhibits  high  optical  sensitivity  and 
high  framing  speed.  It  consists  of  a  photocathode  and 
a  microchannel  plate  (MCP)  in  proximity  focus  with  an 
optional  planar  acceleration  grid  and  an  electro-optic 
crystal  plate,  as  shown  in  Fig.  1.  The  electro-optic  plate 
carries  a  high-resistivity  dielectric  mirror  on  one  side 
and  a  transparent  conducting  electrode  on  the  other. 

In  the  electron-deposition  write  mode,  incoherent  or 
coherent  write  light  incident  upon  the  photocathode 
creates  an  electron  image,  which  is  amplified  by  the 
MCP  and  proximity  focused  onto  the  surface  of  the 
dielectric  mirror.  The  resulting  spatially  varying 
electric  field  modulates  the  refractive  index  of  the 
electro-optic  plate.  Thus  the  readout  light,  which 
makes  a  double  pass  through  the  crystal,  is  spatially 
phase  or  amplitude  modulated,  depending  on  the  crystal 
cut  and  readout  scheme  (polarization  or  interferomet¬ 
ric)  employed. 

The  image  is  erased  by  appropriately  adjusting  the 
device  voltages  and  flooding  the  photocathode  with  light 
so  that  electrons  are  removed  from  the  mirror  surface 
by  secondary,  electron  emission.  Alternatively,  the 
device  can  be  operated  in  the  reverse  mode,  in  which  a 
positive-charge  image  is  written  by  secondary  electron 
emission  and  erased  by  adding  electrons  to  the  crystal. 
For  Pockels-effect  crystals,  this  leads  to  contrast-re¬ 
versed  images.  Further  details  of  the  principles  of 
operation  of  the  device  can  be  found  elsewhere.1 

The  results  described  in  this  Letter  were  achieved 
with  a  vacuum-demountable  MSLM  employing  an 
oblique-cut  LiNbC>3  crystal  plate  330  ^im  thick.  The 
crystal  was  cut  such  that  the  normal  to  its  surface  was 
located  in  the  plane  containing  the  z  and  -y  axes  and 
made  an  angle  9  of  55°  with  the  z  axis.  This  crystal  cut 
optimizes  the  half-wave  voltage  and  the  half-wave 
surface-charge  density  in  the  longitudinal  mode.2 

For  the  general  crystal  cut  with  orientation  9,  it  can 
be  shown  that  the  voltage-induced  phase  retardations 


along  the  induced  dielectric  axes  (x'  =  x  and  y')  in  the 
transmission  mode  are  given  by 

Ai i>x(9,  V)  *  Vi7o3(f*i3  cos  9  +  r-w  sin  9) 

m^Vrio3rx(6)  (1) 

and 

A (0,  V)  =  — Vj7e3(0)[(r13  cos  9  -  roo  sin  9)  cos2  9 

+  (r*33  -  2r42)cos  9  sin2  9] 

•jVi uH9)ry(8),  (2) 

where  Ve(9)  is  the  extraordinary  refractive  index  for 
light  propagation  along  z',  the  normal  to  the  crystal 
surface,  and  polarized  in  the  y'  direction.  It  can  be 
shown  that  ife(9)  is  given  by 

r)e(9)  =  ,  2  .  (3) 

(i?o2  sin2  9  +  i)E-  cos2  9)lu 
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In  Eqs.  (1M3),  V’  is  the  applied  voltage.  <\  is  the  wave¬ 
length  of  the  light  in  vacuuo ,  rj0  and  i)e  are  the  ordinary- 
refractive  index  and  the  principal  value  of  the  extraor¬ 
dinary  refractive  index,  respectively,  and  r13,  r2o,  r33, 
and  r42  are  elements  of  the  electro-optic  tensor  for 
LiNbOs. 

The  relative  phase  retardation  T  is  given  by  T  =  d>y 
-  <t>x,  and  the  corresponding  transmission  half- wave 
voltage  Vtt  is  given  by 

Vxt(Q)  *  — iTm — 73T - 3 — 7m  ’  14) 

TleH°)ry'{0)  ~  r)o3rx(9) 

The  values  of  the  refractive  indices  and  electro-optic 
coefficients  for  LiNb03  place  the  minimum  value  of  V  .  - 
near  J  =  55°. 

It  is  clear  that,  when  this  crystal  is  read  out  in  co¬ 
herent  light  between  crossed  polarizers  (oriented  at  45° 
to  the  x  axis),  the  signal  is  modulated  in  both  amplitude 
and  phase.  Amplitude-only  modulation  can  be 
achieved  by  a  sandwich  of  two  appropriately  oriented 
crystals  of  equal  thickness.  At  6328-A  wavelength,  for 
6  =  55°,  the  calculated  value  of  VxT  is  2300  V,  compared 
with  6200  V  (T  *  <t>x  -  <t>y  -  ir)  for  the  z-cut  material. 
For  crystals  operated  in  the  reflex  mode,  such  as  those 
in  the  MSLM,  VrR  =  Vxt/2,  and  the  half-wave  sur¬ 
face-charge  density  <xTfl  is  given  by  o*r  =  <ifo  VrR/lx, 
where  <o  is  the  permittivity  of  free  space  and  e :  is  the 
dielectric  constant  for  LiNbOs  parallel  to  the  z'  axis. 
The  measured  half-wave  value  of  \\r  for  our  samples 
was  1250  V.  This  corresponds  to  a  calculated  value  of 
of  0.14  MC/cm2. 

The  half-wave  write  time  of  the  standard  MSLM  in 
the  electron-deposition  write  mode  is  given  by 

tiw  *  ^jrfi/J/JU,  (5) 

where  Jpw  is  the  primary-write  current  density  incident 
upon  the  dielectric  mirror.  Note  that  Jpw  will  be  less 
than  the  MCP  output-current  density  Jo  in  devices 
employing  an  acceleration  grid  because  the  grid  inter¬ 
cepts  a  portion  of  J0.  Thus  Jpw  <  J0  « I,/ A,  where  Ia 
and  A  are,  respectively,  the  strip  current  and  the  active 
area  of  the  MCP. 

The  corresponding  half-wave  erase  time  is  given 
by1 

t„  -  2=®  »  £=*  (S'  -  l)"1,  (6) 

Jg  Jpe 

where  the  gap  current  Jg  depends  on  the  effective  col¬ 
lected  secondary  electron-emission  coefficient  S'  and  Jp e 
is  the  primary-current  density  incident  upon  the  di¬ 
electric  mirror.  The  coefficient  S’  depends  on  the  grid 
and  mirror  voltages,  the  energy  distribution  of  the  MCP 
output  electrons,  and  the  electronic  properties  of  the 
dielectric  mirror  materials.  Additionally,  S'  is  often 
much  smaller  than  omax,  the  peak  value  of  the  standard 
secondary  electron-emission  curve,  which  is  measured 
on  an  uncharged  surface  with  a  gap  field  of  0  V/cm. 
Note  that  Jpe  2  Jpw.  The  inequality  holds  when  the 
grid  potential  is  greater  than  the  surface  potential  of  the 
dielectric  mirror  because  the  slightly  repulsive  field, 
which  must  be  maintained  in  the  gap  during  secon¬ 
dary-emission  erasure,  reduces  the  effective  transmit¬ 
tance  of  the  grid. 


A  vacuum-demountable  device  employing  a  330- 
Atm-thick  oblique-cut  crystal  was  fabricated.  This 
device  employed  a  single  dielectric-mirror  coating  on 
the  crystal  and  was  read  out  between  crossed  polarizers. 
The  MCP  used  had  a  strip  current  of  250  *iA  and  a  gain 
of  1.1  X  104  at  1  kV  of  bias.  The  derice  was  studied  with 
and  without  an  acceleration  grid.  The  grid  had  a  spatial 
frequency  of  30  cycles/mm. 

The  strip-current-limited  framing  rate  of  the 
oblique-cut  LiNb03  prototype  device  described  above 
was  30  Hz  at  full  modulation  depth  and  60  Hz  at  50% 
modulation  depth  with  the  250-pA  strip-current  MCP. 
The  corresponding  half-wave  write  (electron-deposi¬ 
tion)  and  erase  times  were  10  and  23  msec,  respec¬ 
tively. 

The  normalized  light-to-light  coherent  modulation 
transfer  function  (MTF)  was  experimentally  deter¬ 
mined  by  using  a  U.S.  Air  Force  resolution  test  chart  as 
the  input  image,  measuring  /max(ri  and  /mini*'),  and 
substituting  into  the  formula 


l*')  /mini*7) 

M I  r  u - 7 - :  :  —  r — TT  • 

m  /maxi*1)  "b  /min(*') 


(7) 


Here  Imax(i>)  and  / mini*'),  respectively,  stand  for  the 
maximum  and  minimum  intensity  values  at  spatial 
frequency  v  in  the  coherent  readout  beam  and  m  is  the 
spatial-modulation  index  of  the  input  write  light. 

The  measured  normalized  light-to-light  coherent 
MTF  of  the  standard  oblique-cut  LiNb03  prototype 
device  was  found  to  have  a  value  of  approximately  1.9 
cvcles/mm  at  50%  contrast  for  the  330-Aim-thick  crystal, 
a  gap  of  500  pm,  and  no  grid.  This  is  in  good  agreement 
with  the  calculated  value  of  2.0  cycles/mm  at  50%  con¬ 
trast  (10  cycles/mm  at  10%  contrast)  assuming  perfect 
transfer  of  the  optical  image  on  the  photocathode  to  an 
electron  image  on  the  dielectric  mirror.  Similar  results 
were  obtained  when  the  grid  was  inserted  and  the  grid 
voltage  was  maintained  above  2  kV. 

By  exploiting  the  secondary-electron-emission 
characteristics  of  the  dielectric  mirror  and  the  charge- 
readout  intensity  characteristic  of  the  crystal,  several 
image-processing  operations  can  be  achieved  by  direct 
processing  in  the  space  domain  (see  Ref.  1 ).  These  in¬ 
clude  contrast  reversal,  contrast  enhancement,  edge 
enhancement,  image  addition  and  subtraction,  and 
analog  and  digital  intensity-level  thresholding  as  well 
as  binary-level  logic  operations,  such  as  AND,  NAND, 
OR,  NOR,  XOR,  and  NXOR. 

Figure  2  shows  some  of  the  results  of  image  processing 
in  the  space  domain  with  the  oblique-cut  LiNb03 
MSLM.  Note  the  significant  improvements  in  optical 
quality  and  spatial  resolution  over  the  earlier  results  of 
Ref.  1  with  a  0.5-mm-thick  z-cut  LiTa03  crystal. 

Figure  2(A)  shows  the  imaging  of  a  1951  U.S.  Air 
Force  resolution  chart.  The  sequence  of  Figs.  2(B)- 
2(D)  illustrates  analog  and  digital  space-domain  in¬ 
tensity  thresholding.  Figure  2(B)  shows  a  four-level 
image.  Figure  2(C)  shows  the  same  image  after 
space-domain  thresholding  in  which  quadrant  1  was 
reduced  to  the  zero-intensity  level  of  quadrant  3,  leaving 
quadrants  2  and  4  unchanged.  Figure  2(D)  is  the  re¬ 
sultant  binary-level  image  after  the  MSLM  with  a 
stored  image  of  Fig.  2(C)  was  illuminated  uniformly 
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Fig.  2.  Some  of  the  image-processing  operations  achieved  through  space-domain  processing  with  the  oblique-cut  LiNbO ; 
MSLM. 


with  light  to  bring  quadrant  4  to  the  level  of  quadrant 
2  and  thresholding  was  again  employed  to  restore 
quadrants  1  and  3  back  to  zero. 

Figures  2(E)  and  2(F)  illustrate  contrast  reversal.  In 
both  cases  the  images  were  written  by  electron  deposi¬ 
tion.  The  crystal  was  first  biased,  by  means  of  secon¬ 
dary  electron  emission,  with  a  uniform  positive-charge 
background  at  either  a  minimum  in  the  transmittance 
characteristic  [Fig.  2(E))  or  at  maximum  transmittance 
[Fig.  2(F)]. 

Figures  2(G)  and  2(H)  illustrate  contrast  enhance¬ 
ment.  In  Fig.  2(G)  the  image  shown  in  Fig.  2(E)  was 
washed  out  with  background  illumination.  Figure  2( H) 
shows  the  contrast-enhanced  result  after  analog 
thresholding.  Figures  2(1)  and  2(J)  illustrate  edge  en¬ 
hancement  of  a  binary-level  image.  Figure  2(1)  results 
from  Fig.  2(e)  when  the  write-image  exposure  corre¬ 
sponds  to  the  full-wave  surface -charge  density  it.,.; 
similarly,  Fig.  2(J)  results  from  Fig.  2(F)  for  the  full- 
wave  input-image  exposure. 

Figures  2(M)-2(R)  illustrate  the  AND,  NAND,  OR, 
NOR,  XOR,  and  NXOR  operations  of  the  horizontal 
and  vertical  bar  images  shown  in  Figs.  2(K)  and  2(L). 

All  the  images  shown  in  Fig.  2  were  achieved  by  direct 
image  processing  in  the  space  domain  using  the  non- 
linear-readout  intensity  characteristic  of  the  crystal  in 
conjunction  with  the  secondary-electron -emission 
characteristics  of  the  dielectric  mirror.  Except  for  Figs. 
2(0)  and  2(P),  all  the  images  were  achieved  in  a  device 
without  an  acceleration  grid. 


Se\  ral  of  the  significant  advantages  of  oblique-cut 
LiNbO:?  over  z-cut  LiNbO;i  and  LiTaO^  for  intensity 
modulation  have  been  demonstrated  in  this  Letter. 
The  polarization-readout  scheme  used  with  the 
oblique-cut  LiNhO:>  crystal  leads  to  a  more  uniform 
optical  field  across  the  crystal  than  does  the  interfero¬ 
metric-readout  technique  that  was  used  earlier  with  the 
z-cut  material.  The  half-wave  voltage  of  the  oblique- 
cut  material  is  reduced  by  a  factor  of  2.4  over  that  of  the 
z-cut  material,  and  this  translates  into  a  corresponding 
increase  in  optical  sensitivity  and  framing  speed. 

However,  for  phase-only  modulation,  the  z-cut  ma¬ 
terial  may  be  preferred  over  the  55° -cut  material  for 
some  applications  involving  very-low-level  unpolarized 
readout  light  since  it  would  not  be  necessary  to  polarize 
the  readout  beam  for  processing  with  the  z-cut  mate¬ 
rial. 
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A  simple  closed-loop  system  for  estimating  optical  phase,  called  an  interference  phase  loop,  is  presented.  In  this 
system  the  output  intensity  from  an  elementary  interferometric  phase  sensor  (e.g.,  Zernike  phase  contrast,  homo¬ 
dyne,  heterodyne,  polarization,  or  shearing  interferometer)  is  detected  and  used  to  drive  a  phase  modulator  in  the 
path  of  the  wavefront  being  measured.  It  is  shown  theoretically  and  experimentally  that  with  large  gain  this  con¬ 
figuration  ignores  amplitude  fluctuations  and  unambiguously  estimates  phase  at  high  speed  over  a  dynamic  range 
of  multiple-*-  radians.  When  self-interference  (e.g..  Zernike  phase  contrast)  is  employed,  monochromatic  light  is 
not  required. 


Real-time  measurement  of  optical  phase  fluctuations 
<t>i(x,y,t )  is  required  for  such  diverse  applications  as 
studying  refractive-index  variations  in  plasmas  or  gases, 
visualizing  transparent  structures,  measuring  the  at¬ 
mospheric  phase  distortion  to  be  compensated  for  in  an 
adaptive  optical  system,  and  extracting  phase  infor¬ 
mation  from  a  coherent  optical  data  processor.  Ele¬ 
mentary  phase  visualization  techniques,  such  as  Zernike 
phase  contrast,  heterodyne,  homodyne,  shearing,  and 
polarization  interferometry,  suffer  from  a  limited  and 
nonlinear  dynamic  range  of  ir  radians  as  well  as  from 
2nir  and  phase-quadrant  ambiguities  in  the  inferred 
phase.  In  addition,  they  do  not  by  themselves  distin¬ 
guish  wavefront-amplitude  variations  from  phase 
variations.  A  variety  of  quite  sophisticated  electrical 
and  optical  systems  has  been  designed  to  overcome 
these  difficulties  and  accurately  estimate  phase.1'5 

This  Letter  presents  a  simple  closed-loop  phase- 
estimation  system  (Fig.  1)  called  an  interference  phase 
loop  (IPL),  in  which  the  output  from  one  of  the  ele¬ 
mentary  interferometers  mentioned  above  is  used  to 
drive  a  phase  modulator  in  the  path  of  the  wavefront 
being  measured.  The  spatially  parallel  drive  signal  to 
the  modulator  is  shown  to  be  the  phase  estimate.  In 
many  cases  the  interferometer  intensity  itself  (e.g.,  at 
C  in  Fig.  1)  also  becomes  an  intensity  representation  of 
the  phase  to  be  measured.  It  is  shown  both 
theoretically  and  experimentally  that,  with  targe  system 
gain,  this  basic  feedback  configuration  can  inherently 
ignore  amplitude  fluctuations,  eliminate  phase-quad- 
rant  ambiguities,  and  continuously  estimate  phase  over 
multiple-T  radians  of  dynamic  range.  When  self-in¬ 
terference  is  employed  (e.g.,  shearing  or  polarization 
interferometry  or  Zernike  phase  contrast),  monochro¬ 
matic  light  is  not  required. 

A  phase-measurement  system  employing  an  ac  in¬ 
terferometer  with  feedback  from  a  phase  modulator  in 
its  local  oscillator  (LO)  beam  was  previously  reported.5 
However,  instead  of  exploiting  the  inherent  properties 
of  optical  phase  feedback,  which  we  present  here,  that 


system  employed  complicated  analog  and  digital  elec¬ 
tronics.  Additionally,  the  IPL  does  not  use  an  LO  beam 
in  the  self-interference  implementations. 

Similar  principles  are  also  implicit  in  adaptive-opti¬ 
cal,  transmitter-predistortion  systems  with  feedback 
derived  from  a  target  reflection  [e.g.,  coherent  optical 
adaptive  techniques  (COAT)1  or  Cathey  et  a/.’’]  and 
closed-ioop  adaptive  receivers.1-7  However,  most  of 
these  systems  employ  heterodyne  or  synchronous  de¬ 
tection  and  more-complex  electronics  to  separate  the 
phase  from  amplitude  variations.  The  IPL  is  also 
somewhat  related  to  a  recently  reported  closed-loop 
bistable  Fabry-Perot.8 

The  homodyne-heterodyne  implementation  will  be 
analyzed;  the  analysis  is  similar  for  self-interference.  In 
Fig.  1  the  modulator  adds  a  constant  phase  shift  0oU.>') 
and  an  electrically  controlled  phase  <t>m(x,y,t)  to  the 
input  wave.  The  modulated  beam  is  assumed  to  suffer 
an  attenuation  a  and  phase  shift  <t>\(x,y)  in  traveling 
from  the  modulator  to  the  detector.  When  the  spatial 
dependence  is  suppressed,  LO  field  Ei  and  modulator 
field  Em  combine  to  produce  the  intensity: 

/  *  /o  +  /i  cos[0,  -  (< i>m  +  <t> o)  -  ( <f>i  +  tt/2) 

-  0i  -  («,  -  w/)f  J,  (la) 

where 


,  \  bE-M 
e,  *  ;  E  ,•  i  e  fc'*  ’ 


"  1 

Ml/.- 

-  1/  .  i.*'  1  7  ♦  : 

A  VP 

P. 

w  -  ♦ 

c,  |L.,h-  . 

t'V'l'  »'  !v  •* 

■ ■«?'.  ,r,'r  <: 

Fig.  1 .  Homodyne-heterodyne  implementation  of  the  IPL. 
MOD  is  a  spatial  phase  modulator.  DET  is  a  detector  array. 
AMP  is  an  amplifier  array.  Em  and  Et  are  defined  at  the 
detector  face. 
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/0  =  «-!£,!-+  Inl¬ 

and 

/i-2a|£,-||£,|. 

This  is  more  conveniently  expressed  as 
I  =  /o  +  h  sin  <t>e. 

Here 


0e  ~  0 i  0m  06 


lib) 

(1c) 


(2a) 


(2b) 


is  the  loop  phase  error,  with  06  =  0o  +  0i  +  (w,  — 
<*>/)f. 

The  detector  in  Fig.  1  can  be  modeled  as  exhibiting 
threshold  and  saturation  by 


0 


io 


\8o(I-It) 
go(/s  -  It) 


I  <  It 
/t  </ < 
I*  I, 


(3) 


The  phase  modulator  is  also  assumed  to  have  a  limited 
dynamic  range;  assuming  that  0  <  0m  <  0S, 


0m 


05 


i  i  <  0 

0  <  l'i  <  0s/g2- 
i  l  -  05/52 


(4) 


The  frequency  response  of  the  detector  and  modulator 
can  be  lumped  into  the  amplifier  transfer  function  H(s); 
thus,  in  the  Laplace  domain, 

Ii(s)  =  gitf  (s)Io(s).  (5) 

Combining  Eqs.  (2)— (5)  results  in  an  implicit  expression 
for  the  modulator  phase  in  its  linear  operating  region: 

$m(s)  =  5o5i52H(s)(/0  -It  +  1 1  sin  4>e).  (6) 

Introducing  the  new  lumped  variables  Go  and  G  t  and 
using  Eq.  (2b), 

$m  =  H(s)[Go  +  Gi  sinf^1;  —  <I>m  —  4*6 )].  (7) 

In  order  to  proceed,  a  single-pole  loop  filter,  H(s)  = 
A/(s  +  A),  will  be  assumed;  infinite-bandwidth  (A  -*  ® ) 
or  integrating  [ H(s )  =  1/s]  filters  can  be  viewed  as 
special  subcases.  Owing  to  system  nonlinearities,  it  is 
easier  to  analyze  the  loop’s  operation  in  the  time  do¬ 
main.  With  the  single-pole  filter  and  Laplace  trans¬ 
form  pair  s  **  d/d t,  Eq.  (7)  results  in  the  equation 

0m  *  A[— 0m  +  Go  —  Gi  sin(0m  ~  0i  +  0i,)J,  (8) 

where  0m  *  d0m/dt. 

Graphical  techniques9  provide  a  powerful  tool  for 
studying  these  transcendental  differential  equations, 
particularly  with  such  complications  as  saturation  and 
threshold.  The  modulator  state-space  plot  (0m  versus 
0m)  of  Fig.  2(a)  is  a  graphical  representation  of  Eq.  (8). 
Static  background  phase  (0s  =  0,  i.e.,  homodyne  or 
self-interference)  is  assumed,  making  it  convenient  to 
reference  the  input  phase  to  0s;  i.e.,  0j  =  0,  -  <t>b  (<t>.  = 
0  when  0;  »  0s).  Limited  modulator  dynamic  range  (0 
5  0ffl  <  0,)  and  detector  threshold  and  saturation  are 
also  assumed.  Whenever  the  derivative  is  positive  (0m 
>  0)  the  modulator  phase  increases,  following  a  trajec¬ 


tory  to  the  right  in  Fig.  2(a);  whenever  0m  <  0  the  tra¬ 
jectory  goes  to  the  left.  Stable  equilibria  (circled 
points)  occur  where  the  trajectories  converge  from  above 
and  below  onto  the  0m  =  0  axis. 

According  to  Eq.  (8),  the  whole  sinusoidal  trajectory 
of  Fig.  2(a)  shifts  laterally  in  proportion  to  0,.  The 
stable  equilibrium  points  of  0m  thus  follow  the  trajec¬ 
tories  in  Fig.  2(b),  continuously  tracking  0,  variations 
within  the  system  bandwidth  (0,  <  A 0;).  These  equi¬ 
librium  trajectories  can  also  be  seen  by  setting  0m  =  0 
in  Eq.  (8);  it  then  follows  that 

0/  =  0m  +  2nir  +  sin-l[(G0  -  0m)/Gi].  (9) 

Thus,  with  large  gain  (Gi  -*  ®),  the  modulator  phase 
estimates  the  input  phase:  0m  *  <t>\  +  2n7r. 

As  0,  increases,  0m  eventually  saturates;  then  0m  = 
0,  and  0m  =  0  until  an  unstable  equilibrium  of  Fig.  2(a) 
reaches  0., .  F rom  there  the  system  follows  the  trajec¬ 
tory,  which  causes  0m  to  jump  to  the  next-lowest  stable 
equilibrium  point.  Similarly,  with  decreasing  phase, 
0m  does  not  stay  pinned  at  its  minimum  value  but 
eventually  returns  to  its  phase-tracking  region  [Fig. 
2(b)]. 

In  practice  there  need  not  be  a  2 nir  degeneracy  in  0m, 
since  the  value  of  the  modulator  phase  before  the  system 
is  turned  on  determines  which  equilibrium  point  in  Fig. 
2(a),  and  hence  which  trajectory  in  Fig.  2(b),  is  followed. 
For  example,  if  a  uniform  calibrated  intensity  is  occa¬ 
sionally  shined  into  the  detector  to  reset  <t>m  at  2tt,  then 
all  the  points  in  the  wavefront  will  be  operating  on  the 
0m  =  0i  +  2tt  trajectory  of  Fig.  2(b). 

Wavefront-amplitude  fluctuations  cause  the  offset 
Go  and  gain  G 1  to  vary.  Variations  in  G0  cause  the 
trajectories  in  Fig.  2(a)  to  shift  up  and  down  between  the 
threshold  and  saturation  limits;  with  large  gain,  only 
small  perturbations  in  the  equilibrium  value  of  4>m  re¬ 
sult.  Provided  the  gain  G\  remains  large,  amplitude- 
induced  fluctuations  in  gain  are  also  inconsequential; 
only  the  error,  which  is  already  very  small,  will  vary. 

Notice  that  the  interferometer  intensity  is  both  an 
unamplified  version  of  the  phase  estimate  [Fig.  1  with 
H(s)  *  1]  and  a  representation  of  0e  [Eq.  (2)].  Thus, 
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Fig.  2.  (a)  Modulator  state-space  trajectory  {<t>m  versus  <im ) 

from  0m  *  M~0m  +  Go  -  G\ sin(0m  -  0,)[  with0,  *  0.  G0  = 
1.8jt,  G 1  =  7.3jt,  and  0,  =  3.8ir.  Stable  equilibria  are  circled, 
(b)  Phase  estimation  performance  (0m  versus  0,). 
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Fig.  3.  Single-spatial-element  homodyne  IPL  implemented 
with  a  Mach-Zehnder  interferometer  and  with  direct  electrical 
readout  of  intensity  |£,(f  )|2  and  phase  <t>i(t). 


Fig.  4.  Experimental  results:  /  is  interferometer  intensity. 
>!>e  is  derived  from  Eq.  (2b).  <bm  is  the  modulator  phase,  i.e., 
the  phase  estimate  of  &.  (a)  Phase  tracking  (G\  *  lOir).  (b) 
Closed-loop  and  open-loop  response  to  wideband  (0-150-Hz) 
phase  fluctuations  with  G  i  =»  20ir  and  two  different  loop-filter 
bandwidths  (Xt  ~  104  sec-1,  X2  *  50  sec-1),  (c)  Immunity  to 
amplitude  |£,|  fluctuations  (Gt  *  20 7r). 


in  analogy  with  a  classical  type-zero  control  system,  the 
error,  sin <t>e,  is  amplified  to  drive  the  output,  4>m.  It  can 
also  be  shown  that  the  input  phase  range  corresponding 
to  the  full  interferometer  intensity  range  of  2 1\  [Eq.  (2)] 
is  extended  from  v  radians  in  the  elementary  interfer¬ 
ometer  to  2Gi  +  it  radians  in  the  IPL  system.  The 
phase  image  (e.g.,  signal  at  C  in  Fig.  1)  could  be  recorded 
with  high-contrast  film  or  a  detector  array  having 
threshold  and  gain. 

This  analysis  has  been  extended  to  include  the  het¬ 
erodyne  implementation,  other  loop  filters  [e.g.,  H(s) 
=  1,  K/S,  or  K'/(s  +  Xi)(s  +  X2)],  and  the  system 
transient  response.  The  IPL  system  is  much  like  an 
electronic  phase-lock  loop,9  particularly  in  the  hetero¬ 
dyne  implementation. 

The  single-element  homodyne  IPL  system  of  Fig.  3 
was  constructed.  A  silicon  photodiode  and  a  transverse 
LiNbO:)  phase  modulator  were  employed.  The  mod¬ 
ulator  had  a  half-wave  voltage  of  200  V  and  a  breakdown 
voltage  >  100  kV;  however,  the  maximum  amplifier 
voltage  (700  V)  limited  the  modulation  to  3.5*-  rad. 
The  system  had  a  dominant  single  pole,  which  was  ad¬ 
justable  from  4  to  4000  Hz.  With  a  high-breakdown- 
voltage  (>400  V)  photodiode,  it  should  be  possible  to 


simplify  the  system  further  by  eliminating  the  ampli¬ 
fier. 

Figure  4(a)  demonstrates  the  phase-tracking  behavior 
of  the  phase  estimate  <f>m  (i.e.,  modulator  voltage)  and 
the  interferometer  intensity  I  (i.e.,  detector  signal  i<j) 
for  large  phase  variations  produced  by  tilting  a  glass 
plate.  The  phase  corresponding  to  a  given  tilt  was 
calibrated  by  opening  the  loop  and  interpreting  the 
nonlinear  output  of  the  resulting  unmodified  Mach- 
Zehnder  interferometer.  Figure  4(b)  displays  IPL 
operation  in  the  presence  of  broadband  (0-150-Hz) 
phase  fluctuations.  Notice  that  with  large  bandwidth 
and  high  gain  all  the  fluctuations  are  tracked  by  4>m  and 
"compensated”  out  of  I.  Compensation  occurs  because, 
with  large  gain,  only  a  small  change  in  I  is  required  to 
drive  a  large  change  in  4>m.  In  Fig.  4(c)  an  electro-optic 
amplitude-only  modulator  varied  |£,|;  note  that  the 
phase  estimate  c pn  is  insensitive  to  intensity  variations 
of  more  than  an  order  of  magnitude. 

Besides  phase  measurement,  the  IPL  system  is  well 
suited  for  other  applications.  For  example,  the  bistable 
(or  multistable)  operating  characteristic  of  Fig.  2(b) 
could  be  exploited  to  construct  a  two-dimensional  op¬ 
tical  digital-processing  element.  In  addition,  Fig.  1  can 
be  viewed  as  an  adaptive  optical  system  to  compensate 
for  atmospherically  induced  phase  distortions,  <t>i(x,y,t). 
With  high  gain  the  wavefront  of  E0  at  A  in  Fig.  1  has 
phase  <bb(x,y),  which  is  a  system  constant  and  can  be 
made  plane  or  any  other  desirable  shape.  The  config¬ 
uration  of  Fig.  1  can  also  generate  the  wavefront  con¬ 
jugate  of  Oi(x,y)  in  a  second  laser  beam  injected  at  A. 
An  IPL  system  with  19  discrete  hard-wired  elements 
and  a  high-resolution  “all-optical”  version,  employing 
a  monolithic  optically  driven  spatial  phase  modulator,1.0 
are  under  construction  and  will  be  presented  in  a  future 
paper. 
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An  “all-optical”  approach  for  the  realization  of  adaptive  optical  systems  potentially  containing  in  excess  of  a  mil¬ 
lion  spatial-resolution  elements  is  reported.  A  phase-measurement  and  -compensation  technique  called  an  inter¬ 
ference  phase  loop  is  employed  in  conjunction  with  a  monolithic  optically  addressed  spatial  light  modulator  (SLM). 
Wave-front  phase  compensation  and  shaping,  and  the  ability  to  ignore  amplitude  fluctuations  and  compensate 
phase  in  real  time  over  multiple  t  radians  of  dynamic  range,  were  demonstrated  in  two  discrete-channel  laboratory 
test  systems  containing  one  and  nineteen  resolution  elements.  All-optical  phase  compensation  with  a  monolithic 
SLM  was  also  successfully  demonstrated. 


High-resolution  adaptive  optical  systems,  containing 
thousands  to  millions  of  spatial-resolution  elements, 
offer  the  potential  of  considerably  improving  perfor¬ 
mance  in  such  applications  as1*3  optical  communica¬ 
tions  through  low-visibility  or  turbulent  atmospheric 
conditions,  imaging  through  the  atmosphere,  and 
high-energy  laser  propagation.  Most  previous2  adap¬ 
tive  optical  systems  consisted  of  arrays  of  individual 
phase  estimators  hard-wired  to  discrete  modulators,  a 
technology  that  cannot  practically  be  extended  to  the 
required  high  resolutions.  Instead,  an  all-optical  par¬ 
allel  processing  approach  is  presented  here. 

Figure  1  is  a  conceptual  diagram  of  an  all-optical 
phase-compensated  communications  receiver.  The 
distorted  received  wave  front  passes  through  an  opti¬ 
cally  addressed  monolithic  spatial  phase  modulator, 
which  adds  the  appropriate  phase  to  present  a  com¬ 
pensated,  approximately  planar  wave  front  to  the 
communications  receiver.  Part  of  the  compensated 
wave  front  is  fed  back  to  a  phase  estimator  whose  output 
is  an  optical  intensity  proportional  to  the  phase  of  the 
compensated  beam.  This  optical  control  signal  drives 
the  spatial  light  modulator. 

It  should  be  stressed  that,  unlike  the  system  pre¬ 
sented  here,  high-resolution  nonlinear  (e.g.,  denerate 
four-wave  mixing)  compensation  schemes4  generally 
require  a  double  pass  through  the  distorting  medium. 
For  example,  a  probe  beam  must  initially  be  sent  from 
the  receiver  to  the  transmitter,  or  a  probe-beam  glint 
must  be  received  off  the  object  in  imaging  applica¬ 
tions. 

Most  existing3  techniques  for  phase  measurement 
appear  to  have  undesirable  performance  limitations  or 
seem  too  complicated  for  implementation  in  the  all- 
optical  phase  estimator  of  Fig.  I.  For  example,  ele¬ 
mentary  interferometers  based  on  schlieren,  hetero¬ 
dyne,  or  homodyne  techniques  are  all-optical.  How¬ 
ever,  without  additional  signal  processing  they  cannot 
distinguish  wave-front  phase  from  amplitude  variations, 
and  they  suffer  from  a  nonlinear  dynamic  range  of  ir 
radians  as  well  as  from  2 nr  and  phase-quadrant  am¬ 
biguities.  However,  in  an  earlier  Letter3  we  showed 


that  closed-loop  configurations  similar  to  Fig.  1  inher¬ 
ently  simplify  phase  measurement.  More  specifically, 
it  was  demonstrated  that  a  particular  technique,  termed 
an  interference  phase  loop  (IPL),  wherein  the  output 
intensity  from  an  elementary  interferometer  is  fed  back 
to  drive  a  phase-only  modulator  in  the  path  of  the  wave 
front  being  measured,  overcomes  the  above-mentioned 
phase-measurement  problems.  Figure  1  becomes  an 
example  of  an  IPL  when  the  phase  estimator  is  an  ele¬ 
mentary  interferometer. 

Figure  2  illustrates  a  specific  “all-optical”  imple¬ 
mentation  of  the  IPL  employing  a  homodyne  interfer¬ 
ometer  and  a  reflection-readout  monolithic  phase 
modulator,  a  microchannel  spatial  light  modulator 
(MSLM).6,7  In  the  MSLM  the  electron  image  from  the 
photocathode  is  amplified  by  the  microchannel  array 
plate  and  proximity  focused  onto  the  dielectric  mirror. 
The  spatially  varying  electric  field  induced  by  the 
charge  distribution  modifies  the  refractive  index  of  the 
electro-optic  crystal  (e.g.,  LiNbOs),  which  phase  mod¬ 
ulates  the  reflected  readout  beam.  A  half-wave  expo¬ 
sure  sensitivity  of  2.2  nJ/cm2  (X  =*  655  nm),  a  framing 
rate  of  30  Hz,  and  a  spatial  resolution  of  10  line  pairs  per 
millimeter  (lp/mm)  (10%  modulation  transfer  function) 
have  been  achieved  with  current  MSLM’s.7  Quan¬ 
tum-limited  sensitivity,  kilohertz  framing  rates,  mul- 
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Fig.  1.  All-optical  adaptive  phase-compensated  communi 
cations  (direct-detection)  receiver. 
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Fig.  2.  An  “all-optical”  high-resolution  1PL. 


The  modulator  phase,  in  Fig.  2,  is  generally  pro¬ 
portional  to  the  write-beam  intensity,  with  some  tem¬ 
poral  response.  For  example,  with  a  single-pole  low- 
pass  time  response,  the  dynamic  behavior  of  the  mod¬ 
ulator  is  described  by3'5 

(dd>m/dt)  +  u0<pm  =  (jJ0[G0  +  Gi  sin (<?,'  -  <?•*)]•  (3) 

Here  G0  =  g(I„  -  It),  Gx  =  glu  g  is  an  electronic  gain 
factor,  and  /,  is  a  detector  threshold  intensity.  Typi¬ 
cally,  -2 x  <  G0  <  2x  and  2x  <  Gi  <  50x.3 

This  nonlinear  system  has  been  shown3-5  to  be  un¬ 
conditionally  stable,  with  equilibria  (d<pm/dt  =  0)  oc- 


tiple  2x  phase  modulation,  and  spatial  resolution  ap¬ 
proaching  20  Ip/mm  are  ultimately  expected.  The 
writing-  and  reading-light  wavelengths  can  range  from 
the  IR  through  the  UV,  depending  on  the  particular 
choice  of  photocathode  and  electro-optic  crystal. 

Theoretical  and  experimental  studies  have  shown3-5 
that  the  IPL  not  only  largely  ignores  amplitude  fluc¬ 
tuations  but  also  has  no  phase-quadrant  ambiguity  and 
can  continuously  estimate  phase  over  multiple  x  radians 
of  dynamic  range.  When  schlieren  interferometry  is 
employed,  monochromatic  light  is  not  required.  In 
operation,  the  modulator  phase  and  the  elementary 
interferometer  output  intensity  become  representations 
of  the  phase  estimate.  With  large  loop  gain,  the  phase 
error  of  the  modulated  beam  (e.g.,  Ec  in  Figs.  1  and  2) 
can  be  made  arbitrarily  small.3-5  Hence  the  IPL  also 
simultaneously  achieves  phase  compensation. 

The  IPL  can  also  be  employed  for  wave-front  conju¬ 
gation3-5  (e.g.,  by  injecting  a  laser  beam  at  C  in  Fig,  2) 
and  high-resolution  bistability/multistability.3-5  In 
many  respects,  the  IPL  can  be  considered  an  optical 
phased-locked  loop  that  directly  tracks  optical  wave- 
front  phase.  It  is  also  somewhat  related  to  bistable 
optical  devices  incorporating  feedback  from  a  detector8-9 
and  to  a  reported  phase-measurement  system  that 
employed  an  ac  interferometer  with  feedback  by  means 
of  a  phase  modulator  in  its  local  oscillator  beam.10 
Similar  principles  are  also  implicit  in  some  closed-loop 
adaptive  receivers2-11  and  some  adaptive  optical, 
transmitter-predistortion  systems  with  feedback  de¬ 
rived  from  a  target  reflection.2-12 

Only  the  most  basic  operational  principles  of  the 
homodyne  interferometer  implementation  of  the  IPL 
are  presented  here.  A  more  detailed  discussion,  which 
includes  heterodyne  and  schlieren  implementations,  can 
be  found  elsewhere.3 

In  Fig.  2,  the  modulator  adds  a  constant  phase  <j>0  ( x . 
y)  and  an  electrically  controlled  phase  <t>m(x,y,  t)  to  the 
input  corrupted  wave.  When  the  spatial  and  temporal 
dependencies  are  suppressed,  the  modulated  and 
local-oscillator  (LO)  beams  combine  at  the  photocath¬ 
ode  to  produce  the  intensity 

I  »  In  +  /(  sin  (1) 
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Fig.  3.  Single-channel  IPL  experimental  results,  (a)  Phase 
compensation  (<j>e  —  constant)  and  simultaneous  phase  esti¬ 
mation  (<t>m  —  G i  =  lOx.  (b)  Closed-loop  and  open-loop 

response  to  0-150-Hz-bandwidth  input  phase  fluctuations 
with  two  different  IPL  loop-filter  bandwidths  (u>!  =  104  sec-1 
and  «2  =  50  sec-1),  G  i  =  20ir.  (c)  Immunity  of  phase  com¬ 
pensation  ( <t>e  =*  constant)  to  amplitude  |  E,  |  fluctuations,  G  i 
=  20x. 


where 

0*  M  0i  —  <t>m  ~  <t>n  ~  0(  5  <t>i'  ~  <t>m ■  (2) 

Here  /„  -  |  £,  | 2  +  |  E,  | 2.  / ,  -  2 1 E,  1 1 E,  \ .  0, '  =  0,  -  <  0,, 
+  0t),  and  0„  +  <t>t  is  a  static  reference  phase  of  the 
system  li.e.,  0/  ■  0  when  0,  ■  <t>„  +  0/). 


Fig.  4.  Nineteen-element  IPI.  experiments.  'at  Nineteen- 
element  discrete-channel  modulator  array,  ibi  Photographs 
of  Fourier  transform-  aberrated  anti  compensated  waves, 
ic).  id)  Response  "i  a  dilfr.u  mm  limited  detector  in  the 
Fourier  plane  to  u  >  -tatt>  and  -di  .ivnamic  phase  fluctuations. 
tO  is  the  /ero-liitht  level  1 
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Fig.  5.  Phase-compensation  results  with  the  “all-optical” 
MSLM/IPL  system,  (a)  Interferogram  of  the  initial  phase 
distortion  (dynamic  range  of  stir  rad),  (b)  Interferogram 
after  IPL  compensation,  (c)  Long-term  exposure  of  (bj  to 
show  residual  error.  The  fine  fringes  in  the  lower  left  of  (a) 
and  (ci  are  due  to  light  that  reflected  from  the  MCP. 


curring  at  <bm  ■  <p/  -  sin_1[(4>m  -  G0)/G  J  +  2nir. 
Thus,  with  large  gain,  G  i  »  |  G0  -  <t>m  | ,  the  modulator 
phase  estimates  the  input  phase;  i.e., 


a*  4>i'  +  2 me. 


The  n  value  generally  remains  constant,  permitting 
continuous  phase  tracking  over  multiple  ir  radians.3-5 
(In  appropriate  applications,  however,  the  system  can 
be  made  to  exhibit  multistability.)  As  long  as  G  i  »  |  <t>m 
-G0 1.  fluctuations  in  \E,  |  and  \Ei  |  have  a  negligible 
effect  on  <t>m  and  <t>e-3,5 

The  phase  of  the  compensated  signal  beam  (e.g.,  Ec 
in  Fig.  2)  is  <j>c  *  4>t  -  <t>m  -  <t>0  -  <t>f  +  <t>u  With  large  G  lf 
and  <bc  —  <t>i-  The  IPL  thus  compensates  <t>c  to 
the  shape  of  the  LO  wave  front,  which  can  be  designed 
to  be  planar.  Note  that  the  intrinsic  modulator  dis¬ 
tortions  0O  are  removed  from  <pc. 

An  all-optical  phase  compensator  similar  to  that 
shown  in  Fig.  2  was  experimentally  demonstrated. 
Since  the  resolution  elements  of  a  homodyne  IPL  are 
spatially  decoupled,  two  discrete-channel  laboratory 
systems  containing  one  and  nineteen  resolution  ele¬ 
ments  were  also  built  to  investigate  IPL  operation  in 
greater  depth. 

Homodyne  IPL  performance  was  demonstrated 
previously  in  a  single-channel  system  that  employed  a 
Si  photodiode  to  drive  a  transverse  LiNbOs  modulator.3 
Figure  3(a)  illustrates  phase  compensation.  The  error 
phase  of  the  compensated  beam,  </>«.  (<t>e  *  <t>i’  ~  <t>m, 
represented  by  the  interferometer  intensity  I)  remains 
relatively  constant  while  <t>m  tracks  >2ir  variation  in  <&■' 
produced  by  tilting  a  glass  plate  in  the  input  beam. 
Figure  3(b)  displays  IPL  operation  with  0-150-Hz- 
bandwidth  phase  fluctuations.  Notice  that  with  large 
bandwidth  and  high  gain  all  the  fluctuations  are  com¬ 
pensated  out  of  <t>i'  {<t>.  at  constant)  and  tracked  by  4>m. 
In  Fig.  3(c),  <t> ,  exhibits  insensitivity  to  intensity  |  £,) 2 
variations  of  more  than  an  order  of  magnitude. 

The  19-element  iest  system  employed  19  Si  photo¬ 
diodes  driving  a  hexagonal  array  of  19  reflective  elec¬ 
trodes  deposited  on  a  disk  of  LiNbOs,  as  pictured  in  Fig. 
4(a).  Images  of  the  Fourier  transforms  (at  the  equiv¬ 
alent  of  plane  A  in  Fig.  1)  of  corrupted  and  IPL-cora- 
pensated  wave  fronts  are  compared  in  Fig.  4(b).  The 
compensated  transform  contains  a  sixfold  symmetric 
pattern  of  repeating  transforms  that  is  due  to  sampling 
of  the  wave  front  by  the  hexagonal  modulator  array. 
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Figures  4(c)  and  4(d),  which  display  the  intensity  seen 
by  a  diffraction-limited  detector  in  the  Fourier  plane 
(proportional  to  the  Strehl  ratio),  illustrate  compen¬ 
sation  (IPL  on)  of  mostly  static  (distorted  Plexiglas  and 
slight  turbulence)  and  dynamic  theated-air)  wave-front 
aberrations,  respectively.  The  compensated  Strehl 
ratio  in  Fig.  4(0  is  increased  by  a  factor  of  7  to  ap¬ 
proximately  the  Strehl  ratio  of  the  unaberrated  wave 
front. 

Figure  5  shows  a  low-resolution  demonstration  of 
phase  compensation  with  an  “all-optical”  MSLM-IPL 
system  similar  to  that  of  Fig.  2.  The  initial  uncom¬ 
pensated  phase  distortion,  which  had  a  dynamic  range 
of  about  tr  radians,  is  depicted  in  Fig.  5(a).  Figures  5(b) 
and  5(c)  display  the  IPL-compensated  interferogram. 
The  long  exposure  of  Fig.  5(c)  reveals  a  very  small  re¬ 
sidual  phase  error  (^ — A/20).  Unfortunately,  the  par¬ 
ticular  prototype  MSLM  that  was  employed  had  defects 
that  prevented  its  use  on  high-resolution  spatial  phase 
aberrations.  Recently  improved  MSLM’s,  however, 
should  permit  truly  high-resolution  MSLM-IPL  phase 
compensation  to  be  realized  in  the  near  future. 

In  conclusion,  a  versatile  wave-front-compensation 
technique  has  been  demonstrated  that  potentially  offers 
the  resolution  as  well  as  the  speed  and  sensitivity  re¬ 
quired  for  a  variety  of  adaptive  phase-compensation 
applications.  The  ultimate  performance  of  this  tech¬ 
nique  will  be  determined  by  the  characteristics  of  op¬ 
tically  addressed  spatial  phase  modulators,  which  are 
continually  improving.  An  MSLM-based  IPL  is  pro¬ 
jected  ultimately  to  be  capable  of  compensating  close 
to  106  (20  lp/mm  over  a  5-cm.aperture)  coherence  cells3 
of  a  highly  corrupted  wave-front  fluctuating  at  kilohertz 
rates. 
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